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PREFACE TO PART II. 

Exemplis melius docemur 
qwam prcBceptis. 

The first part of this work containing problems in 
machine design was pubUshed in 1912, and the favor- 
able reception accorded it gave the author the courage 
to publish the rest of the problems which he had col- 
lected and used in his classes for some years. The 
second part contains problems selected from various 
branches of hydraulics, in the solution of which it is 
necessary to use calculus and analytic geometry. The 
remaining three parts of the work contain problems 
selected from thermodynamics, mechanics of materials, 
and electrical engineering. In this way, a student or 
an engmeer who wishes to review calculus or analytics, 
or to acquire faciUty in applications of higher mathe- 
matics to engineering problems, may select at first the 
part of the work which deals with problems in that 
branch of engineering with which he is most famiUar, 
or in which he is particularly interested. 

The author aimed everywhere to show the student 
not only the possibilities of the mathematical method 
in engineering, but by implication to indicate its limita- 
tions as well. The prominent French engineer Dupuit 
has well said: " Formulae are but tools which may guide 
intelligence, but can never replace it.'' May this serve 
as a warning to both kinds of extremists, those who 
have no use for higher mathematics in engineering 
work and those who follow mathematical deductions 

iii 
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IV PREFACE, 

blindly, forgetting many practical considerations which 
cannot be taken account of in formulae. 

The reader will find the author's views on teaching 
mathematics to engineering students in the preface to 
Part I and in the Dialogue following that preface. He 
is also referred to Part I for a list of reference works on 
mathematics and for an Appendix entitled " What a Sen- 
ior in Engineering ought to know about Mathematics.'' 

The author wishes to acknowledge gratefully the 
assistance of his former student, Mr. Fred G. Switzer, 
M.M.E., who solved most of the problems on hy- 
drauhcs, and of Mr. A. C. Stevens, M.E., instructor in 
Cornell University, who read the proofs. 

Cornell University, Ithaca, N. Y. 
October, 1915. 
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ENGINEERING MATHEMATICS. 



CHAPTER I. 

WATER PRESSXTRE ON SUBMERGED SURFACES. 

A PIPE B (Fig. 1) is connected to a reservoir, A, filled 
with water; the entrance to the pipe is closed with a 
rectangular valve C. The problem is to determine the 
total horizontal force F upon the valve due to the hydro- 
static pressure of the water in the reservoir. 



'. 



m. 



Fig. 1. — A reservoir connected to a pipe. 

The fundamental law which obtains in this case is 
as follows: the force of pressure on an infinitesimal 
area ds at a depth x from the free surface is equal to the 
weight of a vertical coliman of water having ds for its 
base and x for its height. This force does not depend 
on the position of the element ds, being the same for 
vertical, horizontal, or oblique planes, provided they 

3 



Digitized by 



Google 



4 ENGINEERING MATHEMATICS. [Chap. I. 

are all at the same depth x. The plane element, imless 
it is horizontal, must be infinitesimal, as otherwise x 
would be different for different parts of the element. 

This law may be considered to be a combination of 
two simple propositions: (1) The pressure at a point 
in a fluid is the same in all directions; (2) the pressure 
upon a horizontal plane is equal to the weight of the 
colimm of fluid resting upon it. 

Applying this law to our case we find that the force 

acting upon an infinitesimal strip 6 • dc of the valve is 

dF = w 'X 'bdx, (1) 

where w is the weight of a cubic unit of the fluid.* 
Therefore, total pressure 

F = r^wxb dx = ^wb (hi^ - W) . . (2) 
Eq. (2) may be represented in the following form: 



F=^wbih-h)[^^), 



or 

F =W'bk'h, (3) 

where bk is the area of the valve and h is the distance of 
its center of gravity from the free surface of the water. 

Prob. I. — What is the total pressure on a valve 40 cm. wide 
and 30 cm. high, the upper edge being 3 meters below the surface 
of water? Ans. 3750 kg. 

Prob. 2. — A rectangular valve of dimensions h and k (Fig. 1) is 
mounted with its center of gravity h meters below the surface. If 
the height of the water is increased by 10 per cent, what percentage 
change must be made in the dimensions of the valve so that (1) 
the ratio of b and k will not be changed, (2) the total pressure on the 
valve will be the same as before the change in h, 

Ans, 4.66 per cent. 

* For the purposes of this work, consider w constant and for water 
equal to 1000 kg. per cu. m. 
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WATER PRESSURE. 



Expression (3) may be generalized for any shape of 
the valve* To prove this, we note that Eq. (1) holds 
true for any shape of the area if 6 be considered variable, 
so that in general Eq. (2) becomes 

F =w f'xbdx (4) 

But X •bdx represents the static moment of the element 
of area with respect to the free surface of the water. 
Hence, by the definition of the center of gravity, we 
have 

f'xbdx=Sh, . . . • . (5) 

SO that 

F=wSh, .... (6) 

where S is the area of the valve and h is the depth of its 
center of gravity below the free surface. Expressed 
in words, Eq. (6) means that the pressure is equal to the 
weight of a colimm of water having the area of the im- 
mersed surface for its base and the depth of the center 
of gravity for its altitude. 




,1 


't b ». 


;io= 





X 




hr 




t 




k 



Fig. 2. — A submerged trapezoid. Fig. 3. — A rectangle with its upper 

edge at the surface of water. 

Prob. 3. — Find the total water pressure on each side of an im- 
mersed vertical trapezoid (Fig. 2). Consider also the case when 
the lower base is smaller than the upper. 

Ans. F = iwk{2 boh + 2 hih + hohi + bih) . 

Prob. 4. — A vertical rectangle (Fig. 3) has its upper edge at the 
surface of the liquid (^ = 0). At what depth, x, below the surface, 
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6 ENGINEERING MATHEMATICS. [Chap. I. 

should its lower part be cut off horizontally in order to reduce the 
total pressure on the remaining upper part to 1/n of its former 
value? Plot and analyze the curve between x and n. 

Ans. X = hi Vl/n. 
Prob. 5. — Find the value of x in the preceding problem when 
ho is not equal to zero. (This is the general case of Prob. 4. Note 
that when, in the answer to Prob. 5, ho = 0, the answer to Prob. 
4 is obtained.) Hint.—JJoe Eq. (2). 

Ans. X = VV + (^1^ - W)/n. 
Prob. 6. — Solve Prob. 4 when the figure is a triangle with a hori- 
zontal base and its vertex is at the surface. 

Ans. x^hi v^l/n. 
Prob. 7. — Solve Prob. 6 when the base of the triangle is at the 
surface and the vertex is pointed downward. Solution. — Let the 
area of the original triangle (Fig. 4) be S, and let the distance of 

TipgO the unknown division line from the ver- 

^' \ '^ tex be y; the original pressure is F = 

ii V y^"l[ ^^ * i ^i> *^^ *^® pressure on the lower 
, I \/ I triangle is 

Fiq'Z — Asubmerged Fy = wiS* y^W) {h - iy). 

triangle. rpj^^ imposed condition is that F - Fy 

= F/n. Substituting these values of F and Fy, and denoting the 
ratio of y/h by 2, a cubic equation is obtained: 

22»-322+(l-l/n) = (7) 

This equation may be solved by means of Cardan's Rule,* by the 
method of successive approximations if n is given numerically, or 
graphically (see the following problem). 

Prob. 8. — Find graphically the value of 2 in Eq. (7), for n = 5. 
Hint. — Plot t = {Sz^ — 2^) to z as abscissae and find z for which 
/ = 1 - 1/n = 0.8. Ans. z = 0.71. 

Prob. 9. — Find the pressure on a segment of a circle (Fig. 5) 
corresponding to a horizontal chord 2 a and altitude h. Apply the 
answer to the special cases, such as a = 60**, 90®, 120**, and 180**. 
Hint. — Use Eq. (4) in which b and x must be expressed as functions 
of a. Integrate with respect to a. 

Ans. F = ipr* (Ai — r) (a — sin a cos a) + i tor* sin' a. 

Prob. 10. — Same as Prob. 9, except that the segment is part of 
the parabola x = ky^. Ans. F = J to (4 h^hi -3¥). 

* See, for instance, Seaver's Mathematical Handbook, page 23. 
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WATER PRESSURE. 



Formula (6) applies to inclined or horizontal plane 
surfaces as well as to vertical ones (Fig. 6). The pres- 




pe a >A 

Fig. 5. — A submerged segment of a circle. 




Fig. 6. — Water pressure upon an inclined plane area, 
sure on an infinitesimal strip of length b and width dy is 

dF =w 'hdyx, (8) 

so that 

(8a) 



bdy 'X. 



But the expression under the integral sign is the static 
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[Chap. I. 



moment of the area with respect to the free surface of 
the liquid, so that, according to the definition of the 
center of gravity, 

F=^wSh, (9) 

which is identical with Eq. (6). It shows that the total 
pressing on a plane area is the same as long as its center 
of gravity is at the same depth, no matter what posi- 
tion the area itself may occupy with respect to the 
horizontal plane. 

When the position of the center of gravity or the 
pressure F is to be f oimd by integration according to Eq. 
(8a), it is convenient to eliminate x and to integrate in 
the plane AB of the immersed surface. We have 



so that 



X = t/ sm a. 



Jf*Vi 
bydtfj . . 
Vo 



(10) 



where the integral represents the static moment of the 
area with respect to the horizontal water edge at A. 

Prob. II. — A rectangle of area S and height k is hinged on its 

upper side, the hinge being at a 
distance hi below the free surface. 
Its plane is turned into differ- 
ent positions corresponding to the 
limits of the angle a from + 90** 
to —90**. Find the law according 
to which the total pressure F varies 
with a, and sketch the general 
character of the curve F = /(a). 
What is the name of the curve? 
Ans, F = wS{hi + iksma). 

Prob. 12. — Find the total pres- 
sure on a valve made according 
to Fig. 7 and placed at an angle a = 75°, so that the hinged edge 
is 3 meters below the surface of the water. Ana. 12,670 kg. 




Fig. 7. — A circular valve. 
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WATER PRESSURE. 



Center of Pressure. 

In some hydraulic problems it is of importance to 
know not only the value of the total pressure of the 
liquid on an immersed surface, but also the point at 
which the resultant pressure is applied. This point 
is called the center of pressure. A little consideration 
will show that this point always lies below the center 
of gravity of the area. Let AB (Fig. 8) be a submerged 
plane area. The total pressure on the surface is equal 
to the weight of a truncated water-cylinder or prism 




Fig. 8. — Submerged plane area showing center of pressure. 

ACDB, in which ordinates perpendicular to AB are 
equal to the vertical distances of the corresponding 
points of AB from the free surface of the water. Thus, 
AC = Aa', BD = jB6', Gg = Gg\ etc. The volume and 
the weight of this cylinder are equal to those of the 
cylinder AabB with parallel bases, provided that the 
plane ab passes through g, which is the projection of 
the center of gravity G of the given area AB. This is 
because the volumes Cog and gDb are equal, and the 
vaKdity of formula (9) for total pressure F is thus 
confirmed. 
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10 ENGINEERING MATHEMATICS. [Chap. I. 

However, the actual distribution of pressure is gov- 
erned by the cyUnder ACDB, and not AabB; the point 
of application of the pressure, H, corresponds to the 
center of gravity H^ of the actual cylinder, and not to 
the center of gravity (?' of the equivalent cylinder. Since 
the actual cylinder of pressure ACDB is always heavier 
in the part corresponding to the deeper immersion of the 
area AB (right-hand side in the figure), it will be seen 
that H is always lower than G. A more rigid proof of 
this proposition follows from the expression (15) given 
below, for the distance of the center of pressure from the 
surface of the Hquid. 

In order to find this distance yp (Fig. 8) to the center 
of pressure, the general method for finding the position 
of the resultant of any parallel forces is used; namely, 
the moment of the resultant force with respect to some 
axis is equated to the smn of the moments of the com- 
ponent forces with respect to the same axis.* 

Let the horizontal water edge through N be such an 
axis. Then the equation of equiUbrium is 

Fyp^TydF, (11) 

or, substituting the values of F and dF from Eqs. (8) 
and (8a), we have 

yp fbxdy = j'bxydy. 
Remembering that x = t/ sin a, we get 

VvT^dy^ T'hy^dy. . . . (12) 

t/yo *^yo 

* In the following theory it is presupposed that the figure is sym- 
metrical with respect to the axis BN which passes through the center of 
gravity, so that the center of pressure lies in the same vertical plane 
with the center of gravity. If the figure is iuis3mmietrical, two or three 
equations of equilibriimi are necessary. This more general case is not 
considered here, since it is of less practical interest and leads to more 
complicated equations. 
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Chap. I.] WATER PRESSURE. 11 

If the shape of the area is given, b may be expressed as 
a function of j/, so that the integration in Eq. (12) may- 
be performed and yp determined. As an example, let 
it be required to find the position of the center of pres- 
sure in the vertical rectangle (Fig. 1). Here b in Eq. 
(12) is constant, and x = y; therefore 

When ho = Oy t/p = f hi. This result is used in hy- 
drauUc engineering, in the design of supporting walls, 
dams, lock gates, etc., where the pressure of water is 
assumed to be concentrated at a distance of one third 
of the depth from the bottom. 

In expression (12) jby dy is the static moment of 

the figure with respect to the horizontal water edge N 

perpendicular to the plane of the paper; / by^ dy is 

the moment of inertia of the figure with respect to the 
same axis N. For many symmetrical figures it is 
easier to find the moment of inertia with respect to 
the axis passing through the center of gravity of the 
figure than for any other axis; moreover, for most 
figures important in practice, the values of the moment 
of inertia with respect to this latter axis are tabulated 
in various engineering pocketbooks. For these reasons 
it is important to establish a relationship between the 
moment of inertia with respect to an axis passing 
through the center of gravity of a figure and the 
moment of inertia with respect to a parallel axis. 
This is done as follows: Let z denote an ordinate of 
the figure with respect to the horizontal axis in the 
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12 ENGINEERING MATHEMATICS. [Chap. I. 

plane of the figure, passing through the center of grav- 
ity G (Fig. 8). We have then for any point, such as 
MyV = yg — Zy Siud dy = — dz; so that the expression for 
the moment of inertia /jv of the figure with respect to 
the axis N may be written thus : 

I by^dy = - / b(yg-zydz= I b (y^ - z^dz, 

2/0 t/yo t/2i 

or, expanding. 

In = Vo^ r°& dz-2yg f\z dz + f^bz^ dz. 

The three integrals in this equation have the follow- 
ing meaning: the first integral represents the area S of 
the figure, the second integral is equal to zero accord- 
ing to the definition of the center of gravity; the third 
one, by definition, is the required moment of inertia 
Iq with respect to the horizontal axis passing through 
the center of gravity. Thus, we have 

lN = SyJ^ + Ia (14) 

In words: the moment of inertia /jv of a figure with 
respect to any axis N in its plane is equal to the moment 
of inertia Iq with respect to the parallel axis passing 
through the center of gravity pliLS the area of the figure 
multipUed by the square of the distance between the two 
axes. Combining now Eqs. (12) and (14) we obtain: 

ypS-ya = SyJ^ + lGy 
or 



fo 



This equation is often convenient for the determina- 
tion of the center of pressure. Incidentally it shows 
that the center of pressure is always below the center 
of gravity, because the last term on the right-hand side 
of the equation is essentially positive. 
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Chap. I.] WATER PRESSURE. 13 

Prob. 13. — Show how to obtain result (13) from Eq. (15). 
Prob. 14. — Find the position of the center of pressure for an 
immersed circle. Solution, — 

7g = 2 Cbz^dz =^ 4 fV Vr^^IT^dz. 

Take angle B (Fig. 7) for the independent variable; then 2 = r cos ^; 
dz = —r sin B dB; hence 

Ig = 4: P Vcos2^.rsme.rsin^d^ = ^ C 8m^2B*d{2B); 
or 

Hence Eq. (15) gives 



lG=i^. 



or 



When the upper edge of the circle touches the free surface of the 
water, yg = r, and 

In other words, the center of pressure lies in this case at a distance 
of i r below the center of the circle. When Vg > r this distance is 
always less than J r. At an infinite depth of immersion the center 
of pressure coincides with the geometric center of the circle. 

Prob. 15. — Apply Eq. (15) to a trapezoid (Fig. 2) with hori- 
zontal parallel sides, when the upper edge is in the surface of the 
water (Ao = 0). Simplify the result for the case of (a) a rectangle, 
(6) a triangle with the vertex up and (c) a triangle with the vertex 
down. 

Prob. 16. — Solve Prob. 15 when h^ is not equal to zero. 
Prob. 17. — Find the center of pressure for the segment in Prob. 9. 
Prob. 18. — Find the center of pressure for the segment in Prob. 
10. 
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CHAPTER II. 
DEPTH OF IMMERSION OF FLOATING BODIES. 

The fundamental law which determines the depth of 
immersion of a floating body is as follows: the weight 
of the displaced Uquid is equal to that of the floating 
body. This is the necessary and evident result of the 
condition that the force of gravity be balanced by the 
upward pressure of the liquid. Let it be required, for 
instance, to calculate the depth of immersion x, of a 
vertical log having a cross section of 0.08 sq. m., 2 m. 
long, and weighing 40 kg. per m. The weight of the dis- 
placed water is x • 0.08 • 1000 kg. and this must be equal 
to the weight of the log, which is 40 • 2 =80 kg. Or, 
X . 0.08 . 1000 = 80, from which a: = 1 m. 

In order to determine the depth of immersion, when 
the floating body has a variable horizontal cross section, 
one has to know the expression for the volmne of the 
immersed part as a fimction of the depth of immer- 
sion. The following problems illustrate the general 
method, and afford practice in calculus and in analytic 
geometry. 

Prob. I. — A solid cone of altitude h floats with ita apex up- 
ivard. The weight of the cone is TF. Determine the depth of im- 
mersion and plot a curve giving the per cent depth as a fimction 
of the specific weight of the material. Solution. — Let x be the 
•distance from the apex of the cone to the surface of the water, and 
let V be the total volume of the cone. This volume is later elimi- 
nated from the solution. The unknown volume of the part above 
the water is v {x/Kf^ because the volumes of similar cones are as the 
cubes of their altitudes. Therefore, the volume of the immersed 

14 
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paxt iav — V {x/hy, and we have the relation wvll — (x/h)^ = TF, 
where w is the weight of a cubic unit of water. Solving for x gives 

x^h v^l - W/iwv). 

The required depth of immersion iah — x and the specific weight 
a of the material of the cone is equal to TF/v, so that the equation of 
the required curve is 

h-x 



V w 



Values of a are used as abscissse, and per cent immersion, {h — x)/h, 
as ordinates. 

Prob. I a. — Solve Prob. 1, for a cone floating with its apex down- 
ward. Am, (h — x)/A = \^S/w. 

Prob. 2. — A prismatic body of specific weight a and having its 
cross section shaped according to an ordinary parabola (Fig. 9) 




Fig. 9. — A floating parabolic cylinder. 

floats horizontally. Prove that the depth of immersion y is deter- 
mined from the equation 

where h is the total height of the body. 

Prob. 3. — Solve Prob. 2 when the cross section of the floating 
body is an ellipse, of which the vertical diameter is 2 A and the hori- 
zontal diameter is 2 6. Deduce an expression for the depth of im- 
mersion. 

Am, irSh/w=l{y-h)/h] V2hy-y^ +h8m''^l(y-h)/h]-\-iirh. 

Prob. 4. — A circular log of diameter 30 cm. and of specific 
weight equal to 0.6 that of water floats horizontally. Find its 
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depth of immersion. Solution, — The area of the segment corre- 
sponding to an angle 2 e (Fig. 10) is r^ (^ — } sin 2 e), so that the 
equation of equilibrium is 1000 r^ (^ - i sm 2 e) = 0.6 X 1000 X ttt^ 
or 

^- isin2^ = 0.6ir (1) 

This equation can be solved only by trial or by plotting a curve 
showing values of [e -- i sin 2 0] and taking the abscissa for which 
the ordinate is numerically equal to 0.6 ». While the angle of im- 
mersion e does not depend upon the radius of the log, the depth of 
immersion x evidently does, because 

x = r(l-cose) (2) 

It is left to the student to solve Eq. (1) and to calculate x from 
Eq. (2). Am, = 99°; x = 17.3 cm. 

Prob. 5, — A hollow cylindrical pontoon of radius r and length 
I weighs p kg. Deduce the equation of the curve showing the rela- 





FiG. 10. — Cross section of a 
floating circular log. 



Fig. 11. — A floating sphere. 



tion between the depth of immersion and the useful load which the 
pontoon can carry. Sketch the general character of the curve. 
Solution, — According to the formula used in the preceding problem, 
the weight of the displaced liquid is wlr^ (^ — i sin 2 0), and this must 
be equal to p + P, where P is the useful load on the pontoon. Angle 
may be eliminated and expressed through the depth of immersion 
X by means of formula (2) of the previous problem. 

Prob. 6. — A sphere of radius r and weight Q is floating in a 
fluid of specific weight w. Find the depth of immersion. Solution. 
— The volume of an infinitesimal horizontal slice (Fig. 11) of radius 
yladV == iry^dxy where y^ = r^ — {r — xY = 2rx — x\ Substitut- 
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ing this value of y^ and integrating the expression for dVto find the 
volume of the inunersed segment, we obtain 

7 = J*^ {2rx - a;2) (to = irx'fr - |\ 

so that the equation of equilibrium becomes 

wrx^lr "" q) ~ Q- 

This is a cubic equation which may be solved by Cardan's rule; in 
numerical appUcations it may also be solved by successive approx- 
imations, or graphically, by plotting expressions si^ and J a;' to a; 
as abscissae. 

Prob. 7. — Solve the preceding problem when the floating body 
is an ellipsoid with semi-axes a, b and c, the axis c being vertical. 
Solution, — The equation of the ellipsoid is 

^^t^t^l (3) 

All horizontal cross sections of the ellipsoid are ellipses. Let 
the horizontal cross section at a distance z from the center have 
semi-axes p and g. Then its equation is 

— 4-^ = 1 (4) 

and the volume of an infinitesimal horizontal slice of thickness dz 
is dt; = irpq dz. The volume of the immersed part is 

F = ^ (""^qdz (5) 

•^ c 

The limits of integration correspond to the case when more than 
one half of the body is under water. If the value of 2 in Eq. (7), 
below, should come out negative, it would show that less than one 
half of the ellipsoid is immersed. 

Before integrating Eq. (5) it is necessary to express p and q 
through 2. Assimiing z in Eq. (3) to be constant, gives the cross 
section corresponding to Eq. (4), so that Eqs. (3) and (4) become 
identical. But Eq. (3) may be represented in the form 



so that 



02(c2-.22) ' 52(c2_22) 

c c 
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Substituting these values in Eq. (5) gives 



-c 

or 



The equation of equilibrium is 

Q = wV, (7) 

from which z can be determined. 

Prob. 8. — Plot a ciure showing the relation between the spe- 
cific weight and the depth of flotation of the ellipsoid considered in 
the preceding problem. Take the dimensions of the ellipsoid to be 
a = 3 m.; 6 = 2.2 m.; c = 1.8 m.; the limits of specific weight to 
be 200 and 1000 kgs. per cu. m., and the Uquid to be water. 

Am. d = 200; 2 = - 0.66. 
5 = 600; 2 = + 0.25. 
5 = 800; 2 = + 0.77. 

Prob. 9. — Show that when a = b = c^ the solution indicated 
in Prob. 7 gives the same result as in Prob. 6. Hint, — Note that 
the limits of integration are taken differently in the two problems. 

Prob. 10. — Solve Prob. 5 for a cylindrical pontoon with hemi- 
spherical ends. , 

Ans, p + P = w{lr^ cos"^[(r — x)/r] — l(r — x) V2rx — x^ 

+ ^2(r-Jx)}. 

Prob. II. — A hypothetical highly compressible fluid has specific 
weight equal to 0.5 on the surface and to 1.6 at a depth of 3 m., 
the density increasing uniformly with the depth. Determine the 
depth to which a parallelopiped of height 1 m. would sink in such 
a fluid, its specific weight being 1.15. Solution, — The specific 
weight at a depth x is 5 = 0.5 + J (1.6 - 0.5) x = 0.5 + 0.367 x. 
Therefore, if the bottom of the parallelopiped is at a depth z it 
displaces a weight of the liquid equal to 



x 



«-i 



Adx^S, 



where A is the horizontal area of the body. This weight must be 
equal to the weight 1.15 A • 1 of the parallelopiped itself. Thus, 



/. 



(0.5 + 0.367 x)d^ = 1.15. 
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After integration, this expression is reduced to 0.5 • 1 + § • 0.367 
(22-1) = 1.15, from which z = 2.27 m. 

Note. — Since the density of the fluid varies in this case accord- 
ing to the straight-line law, the integration, strictly speaking, is 
not necessary, and the law of averages may be directly applied. 
Namely, the average density of the displaced fluid is at the center 
level of the parallelopiped, and this average density must be equal 
to the density of the body itself. Thus, we get 0.5 + 0.367 x = 1.15, 
whence x = 1.77 m. The depth of immersion of the bottom of the 
parallelopiped is therefore 2 = x + 0.5 = 2.27 m., or the same result 
as before. When the density varies according to a more general law 
than the straight-line law, or when the body is of a variable cross 
section not symmetrical with respect to the middle horizontal plane, 
an integration is necessary in order to obtain an expression for the 
weight of the displaced liquid. 

Prob. 12. — A very light sphere of radius r must be immersed 
in the liquid given in the preceding example. Plot a ciu^e show- 
ing the weight W of shot to be put into the sphere for various 
depths of submersion, including also partial immersion. 

Ans. From x = Otox = 2r,W = irx'^ {0.5r + x (0.245r - 

0.167) -0.092x2} ; 
Froma; = 2rtox = oo,W = fxr' {0.5 + 0.367 (x-r)}. 




Fig. 12. — A submerged body of irregular shape. 

Prob. 13. — A given irregular body (Fig. 12) is determined by 
the relation q = f (x), where g is a horizontal cross section and x 
is the distance of this cross section from the bottom of the body. 
The density S of the liquid is given as a function of the depth z 
thus: 5 = (2). It is required to establish a relationship between 
the imiform density D of the body and its depth of submersion. 
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Solution. — The mass of the fluid displaced by an infinitesimal slice 
of the body at the depth zisqS*dz; the mass of the corresponding 
part of the soHd is qD dx. Hence, the equation of equilibrium is 



J*H /•h 

Sqdz = D i qdx. 



'i/-A 



In the left-hand member of this equation q must be expressed 
as a function oi z = H — Xy in order to be able to integrate. Then ■ 
the foregoing equation, becomes 

f" 4> {z)f{H -z)dz = D fV W dx. 

In integrating this equation the depth ff is to be considered con- 
stant; the result is an equation from which the required depth H 
may be determined. 

Prob. 14. — Apply the solution given in Prob. 13 to the case 
when the immersed solid is a body generated by revolving an arc 
of the cubic parabola x = ay^ about the vertical a>-axis, and d = 
(2) = dQ + hz + cz'^. 

Am. 3 chyil - 3 A (6 + 2 c)/8 + 3/5 + 3 8o/5 + 3 Hb/5 

When the floating body is of an irregular shape which 
cannot be expressed by an equation, the volume of its 
immersed part can be calculated only approximately, 
by subdividing it into a large number of thin hori- 
zontal slices. An important case of this kind is that 
of a ship, when it is necessary to calculate its displace- 
ment. Let the hull be subdivided by equidistant hori- 
zontal planes at a distance h from each other, beginning 
from the keel. Let Aq, Ai, . . . , A„ be the areas of 
these cross sections, or water planes, as they are called. 
For the sake of convenience let these areas be plotted 
as ordinates to the distances ft, 2 ft, 3 ft, etc., as abscissae 
(Fig. 13). Connecting the ends of these ordinates by 
a smooth curve one can obtain graphically the areas of 
intermediate water planes. The volume of water dis- 
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placed by the ship at a draft OC is equal to the area 
OABC, because this area represents the expression 



x=0 



A dx. 



In this expression A dx is the volume of an infinitesimal 
sUce of height dx. The integration may be performed 
graphically by means of a planimeter. 




Fig. 13. 



Vertical Distances ^ 

- Areas of water planes plotted against distances from 
the bottom of the ship. 



Multiplying the volume of the displaced water by 
the weight of one cubic unit of water, one obtains the 
weight of the displaced water and consequently the 
weight of the ship itself with its cargo. Repeating 
similar calculations for different values of draft OC, 
a new curve may be plotted giving displacement in 
metric tons to draft in meters as abscissse. From this 
curve it is possible to determine the draft of the ship 
with a given cargo. 

If it is preferred to determine displacements analyti- 
cally, without plotting the curve shown in Fig. 13, a 
certain law must be assumed according to which the 
areas vary between each two consecutive measured 
ordinates. 
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The simplest law is that of a straight line, so that 
each strip in Fig. 13 becomes a trapezoid. The areas 
of these trapezoids are ^ h (Ao + -Ai), § h (Ai + A2), 
. . . , ^ ft (A„_i + An), and their sum is 

F = ft{iAo+Ai+ . . . +A,_i + iAnl. . (8) 

This formula is known as the trapezoidal formula, and 
is accurate enough for many practical applications. 
However, when extreme accuracy is required, or when 

the cross section of the 
immersed body changes 
rather rapidly with the 
depth, areas must be 
computed for a large 
\ number of water planes. 

\ A more accurate formula 
\ which requires fewer 
^^ terms is obtained by 
Fig. 14.— A parabola used to deduce regarding the curve be- 
unpsons e. tween each three consec- 

utive ordinates in Fig. 13 as the arc of a parabola the 
axis of which is parallel to the axis of ordinates (Fig. 
14). It may be proved (see Prob. 17 below) that the 
area of two adjacent strips of equal width h is 

t; = H(a+46+c) (9) 

Therefore, assiuning the number of strips in Fig. 13 
to be even, and adding the above expressions for every 
two consecutive strips, we obtain: 

7 = i A [Ao + 4 (Ai + Aa + A5 + • • • + A„_i) 

+ 2 (A2 + A4 + • • • + A„_2) + AJ. . (10) 

This formula is known as Simpson's rule, and is widely 
used in engineering practice. 



^JV^-^rh- 
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Prob. 15. — The areas of water planes of a ship are as follows: 

Water plane No. 0J1234567 
Areamsq. m. 20 237 333 450 516 556 583 605 620 

Cross sections are drawn every meter; water plane No. corre- 
sponding to the keel, and No. 7 to a draft of 7 m. Calculate the 
volume of displaced water by means of Simpson's formula and de- 
termine the per cent error if formula (8) is used. Note. — Water 
plane No. } is introduced because the area changes more than ten 
times between Nos. and 1. In using Simpson's formula apply 
it between Nos. 1 and 7, and add the area between Nos. and 1 
determined according to formula (9). The same remark applies 
to the use of formula (8). 

Ans. 3415 cu. m.; 0.64 per cent. 

Prob. 16. — Using the data of the preceding problem determine 
the weight of the ship plus its cargo for drafts of 7, 5, 3, and 1 m., 
and plot these weights in metric tons to drafts as abscissae. The 
weight of salt water is 1026 kg. per cu. m. One metric ton is equal 
to 1000 kg. Ans, 3504; 2265; 1129; 222. 

Prob. 17. — Prove formula (9). Solution, — The parabola in 
Fig. 14 is turned into its usual position (Fig. 15) corresponding to 
the standard formula 2/^ = 2 px. The area of an infinitesimal hori- 
zontal strip is (m — x) dy = (w — y^/2 p) dy, so that the area be- 
tween h and c is 

By analogy, for the strip between a and h we have 

Dp 

Adding these two expressions, we get for the total area between a 
and c: 

3 V 2p / 

or, since J {k +f) = g, 
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But the extremities of the ordinates k, g, and /lie on the parabola, 
so that 

Jc^ = 2p(m-c); g^ = 2p(m-b); /* = 2p(w-a). 




Fig. 15. — A parabola used to prove formula (9). 

Substituting these values in the expression for v, we get 

t; = J ^ [6 m — (m — c) — 4 (w — 6) — (w — a)], 

and after reduction we obtain Eq. (9). 

Prob. i8. — Prove Eq. (9) using the equation of the parabola in 
its original position, as in Fig. 14. 
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CHAPTER III. 

TIME OF DISCHARGE OF LIQUIDS FROM 
VESSELS. 

In some practical problems it is necessary to calcu- 
late the time in which a given vessel or a storage reser- 
voir (Fig. 16) may be emptied, or the level of the Uquid 
lo^s^ered to a given height. 

The case of discharge tak- 
ing place through a small 
orifice in the bottom of the 
vessel is theoretically the 
simplest, because (1) the 
pressure of the liquid, or '^ J||, 

the head, is the same at Fig. 16. — A vessel being emptied 
all parts of the opening ; (2) through an orifice in the bottom, 
the orifice being small, the remainder of the fluid may be 
considered at rest. The fundamental law of discharge 
in this case, the so-called Torricelli's theorem, is ex- 
pressed by the equation 

v = V2^, (1) 

where v is the velocity of discharge, h the distance to the 
free surface of the liquid (head), and g is the accelera- 
tion due to gravity. If h is in cm., and g in cm. per 
(second)^, v is in cm. per second. The identity of the 
physical dimensions on both sides of the equation is 
satisfied because 



sec. ^ (sec.)' 

25 



cm. ^^ 

X cm. 
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Formula (1) follows directly from the law of the con- 
servation of energy. When a particle of liquid of mass 
m is at the smf ace of the Uquid in the tank, this particle 
possesses a potential energy mg • h, with respect to the 
bottom of the vessel. During the discharge all this 
energy is converted into the kinetic energy J mv^. 
Equating the two gives 

from which Eq. (1) follows directly. 

Velocity v multipUed by the cross section A of the 
orifice gives the theoretical discharge vA, or the quan- 
tity of hquid passing through the orifice in unit time. 
In reality, the discharge may be considerably smaller, 
due to liquid friction and to a contraction in the issu- 
ing vein. Therefore, the product vA is generally multi- 
phed by an empirical coefficient c, less than imity, and 
the actual discharge is represented by the semi-empiri- 
cal formula 

Q = cAv = cA V2gh (2) 

Here c is the so-called coefficient of discharge, which 
has been determined experimentally for various cases 
of discharge. In all the numerical problems below we 
shall assmne c = 0.60, though in actual practice the 
student should consult an engineering handbook for 
the value of c which applies to his particular problem. 
There are two reasons why the actual discharge is 
less than the theoretical: (a) The issuing particles do 
not move in parallel paths, but form the so-called con- 
tracted vein, so that the actual cross section is less than 
that of the orifice; (&) part of the potential energy 
mgh is wasted in eddies formed near the orifice, so that 
the actual velocity of discharge is less than that given 
by Eq. (1). 
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The simplest case to be considered here is that of a 
prismatic or cylindrical vessel. Let the original level 
of the hquid be h above the orifice; let the final level be 
h2, and let the level be x above the orifice at a time L 
During an infinitesimal interval of time, dt, a quantitjr 
Q dt of the hquid flows out, and as a result the level 
sinks by dx. Hence, if the area of the horizontal cross sec- 
tion of the vessel is F, we have the differential equation: 
cAV2gxdt= -Fdx (3) 

It is important to note the minus sign on the right- 
hand side of the equation. The height x decreases 
with the time, hence dx is a negative quantity. But 
the left-hand side is essentially positive. Therefore, in 
order to make the right-hand side positive, —dx must 
be taken. Separating the variables and integrating, 
we get 

^0 Jhi y/x 

or 

cAV2g't==2F{Vhi-Vh2). . . (4) 

This equation gives the relationship between the quan- 
tities which enter into the problem, and from this 
equation one of the quantities may be determined if the 
rest are known. It will be seen that without the minus 
sign in Eq. (3) the result (4) would have no meaning, 
since hi > /12, and the right-hand side would be negative. 
The time of emptying the vessel is obtained by 
letting /12 = 0. We have then* 

t=-^^ (5) 

cA V 2 ghi 

* It is convenient to preserve the expression V2gh in hydraulic 
formulas, whenever possible, for two reasons: (1) it has a physical 
meaning, being a velocity, while V2g_has none; (2) tables have 
been published which give values of V2 gh for numerical values of A. 
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But cA V2 ghi is the original rate of discharge, and Fhi 
is the capacity of the vessel. Hence, it actually takes 
twice as long to empty a vessel as it would to discharge 
the same amount of .liquid from the orifice were the 
original head to be maintained. In oth^r words the 
average velocity of discharge in emptying the vessel is 
equal to one half the initial velocity. 

Prob. I. — A cylindrical tank, 1.25 m. in diameter, is being 
emptied through a romid opening m the bottom, 5 cm. in diameter. 
The original level of the liquid is 2.5 m. above the bottom. To 
levels of water in the tank as abscissae plot as ordinates the time 
which it takes to reach these levels, and also the corresponding 
instantaneous rates of ^discharge in liters per second. Assume 
c = 0.60. 

Ans, h = 2.5 m., q = 8.2 l./sec, t = sees. 

h = 1.25 m., q = 5.84 l./sec, t = 217 sees. 

A2 = m., g = l./sec, t = 743 sees. 
Prob. 2. — Either curve required in Prob. 1 is a parabola, so that 
instead of calculating its ordinates point by point, it may be easier 
to find the vertex of the parabola and the parameter, and thus 
reduce it to the standard form y^ = 2 px. Show how this may be 
done, and check the results with the curves already obtained. 

One of the important appUcations of the foregoing 
theory is the calculation of the time necessary for empty- 
ing a canal lock. When a ship canal is constructed 
through a hilly country, it is impossible to follow the 
profile of the route, on account of the resulting rapid 
flow of water. Therefore the canal is subdivided into 
sections separated by locks (Fig. 17). A lock consists 
of a chamber sufficiently large to accommodate a ship 
or a barge, and is provided with two gates A and B, 
with sluices or openings C and D. When it is desired 
to conduct a ship from the upper to the lower section 
of the canal, the ship being in position 1, sluice D is 
opened and the lock filled to the level of the upper sec- 
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tion of the canal. When the levels are equaUzed gate 
A is opened and the ship is admitted into the lock. 
Then gate A is closed and sluice C opened to empty 
the lock to the level of the lower section of the canal. 
When this is done gate B is opened and the ship is hauled 




Fig. 17. — A canal lock. 

out of the lock. In transferring the ship from the lower 
to the upper section of the canal the same operations 
are repeated in the reverse order. The time which it 
takes to transfer a ship is consumed mostly in filling 
and in emptying the lock. Therefore it is of great prac- 
tical importance to know how to compute this time, in 
order to determine the capacity of the canal for traffic. 

In applying the foregoing theory of discharge to the 
case of orifice C one must keep in mind that orifice C is 
a submerged opening, so that the hydrauUc head acting 
upon particles of water in it is equal to the difference 
of the heads in the lock and in the lower section of the 
canal. It will be seen that this difference is equal to 
the difference of levels of water and does not depend 
upon the position of the orifice. For this reason, the 
theory is applicable to orifices of any size, and not only 
to small orifices, as in the case of discharge from the 
side of a vessel into the air. 

Thus, formula (5) is directly applicable, provided 
that hi is understood to be the original difference of 
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levels of water, and provided that the lock is prismatic, 
that is, its walls are vertical. If the side walls are 
slanted, cross section F in Eq. (3) is a function of x, and 
this is to be taken into account in performing the inte- 
gration. A lock with a variable cross section is con- 
sidered in the following problem. 

Prob. 3. — Find the time necessary to empty a lock chamber 
with side walls slanting at an angle 6 from the vertical. The length 
of the lock is Z; its width, when full, is 6. The sluice is submerged 
and has an area A, The original difference of levels is h. 

Solution. — Let the difference of levels at a time t be x, and let 
bg equal the corresponding width of the free surface of water. 
Then according to formula (3), 

cA V2 gx dt = —Ibx dx, 

where, according to the geometric shape of the lock, 

6, = 6 - 2 (A - x) tan d. 

Separatiag the variables and integrating, we get 

cAV2g.t = - Z r[6 - 2 ih-x)tme]/y/x*dx, 

or 

cAV2g*t = 2lVh\b-2h+ihtQjie], 
whence 

^ 2lh\b-2h + ihiand\ .^. 

cAV2gh 

Prob. 4. — When the slant of the walls of a lock is small, or the 
walls are nearly vertical, it may be sufficiently accurate to use 
formula (5) instead of (6), taking for F the average area of the free 
surface of water between its highest and lowest positions. Deduce 
an expression for per cent error, as a function of the ratio b/h and 
of the angle d, when formula (5) is used, instead of (6). 

Ans. 100X ,,,,^-j^,; , . 
(b/h) -2 + itajie 

Prob. 5. — Represent per cent error deduced in the preceding 

problem by a series of curves for values of b/h from 20 to 2, and for 

values of batter from 1 : 20 to 1:1. Analyze and discuss the 

curves. Note, — A batter of 1 : 20 means that the surface of the 
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wall rises 20 feet for every foot of horizontal projection; in other 
words, tan $ = 0.05. 

Prob. 6. — A vertical cylindrical vessel of horizontal cross sec- 
tion F is filled with a Uquid to a level hi. Two small holes of area m 
each are pierced simultaneously in its side, one at a height a, the 
other at a height 6, above the bottom. Determine the time in 
which the level of the liquid will sink to a level M- Consider two 
cases: (A) ^ > 6 > a, and (B) b> h2> a. Check the answer 
by applying it to some simple case, for example a = b. 

Arts. (A) t 2^!_-_,{(A,-a)*-(^-a)^-(Ai-6)* 

3(b-a)cmV2g . 

+ (^ - oyi ; 

(B) t = ——^^^-=Hhi^a)i+2(]b^a)^^(h^b)^ 
3(b-a)cmV2g , 

— 3(6— a)(A2— a)*}. 

Prob. 7. — A swimming pool is 20 m. long and 9 m. wide. The 
depth of water varies uniformly from 100 cm. to 250 cm. What 
size opening must be made in the bottom of the deep end in order 
to empty the pool in 15 minutes? Hint, — For the upper 100 cm. 
use formula (4) ; for the rest integrate Eq. (3) in which F is variable 
and is determined from the proportion F -^ (20 X 9) = x 4- 150. 

Am. 1150 sq. cm. 

Prob. 8. — Show that when water is discharged from a hollow 
pyramid or a cone with the apex down, through a small orifice in 
the apex, the time of emptying is 20 per cent longer than the time 
which would be required to discharge the same volume at the 
initial head. Hint. — In formula (3) the variable cross section is 
determined from the proportion Fg -r- Fi = x^ -t- h^, 

Prob. 8a. — A water clock is an instrument which consists of 
a vessel with a small hole in the bottom. The liquid which fills the 
vessel is allowed to escape through this hole, and time is measured 
by the position of the free surface of the liquid. Show that, in order 
to have uniform divisions, the vessel must be so proportioned that 
areas of the surface of the liquid vary as square roots of the vertical 
distances to the orifice. 

Prob. 9. — A reservoir has such a shape that the area F^ of 
horizontal cross sections varies'according to the law Fx = A. + Bx + 
Cx^, Find the amount of time nec^sary to lower the level from 
^1 to ^ through an orifice of area A made in the bottom. 
Ans. t^F[2A{Vhi- v%5)+jB(M-^')+«C(;ii*-^*)] 

/{cA V2^). 
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Prob. 10. — A spherical bowl of radius r is provided with a small 
opening of area A at the lowest point. Plot a curve showing the 
time necessary for emptying the bowl, as a function of the original 
level, from a; = 2 r to a; = 0. Take the time necessary to empty 

the semi-sphere -as a unit. 

Ans. t = IT (frx* — ix*)/(4cA ^2g); when a; = r, f = 1, 

so that A = 15/(14 r*). 

Prob. II. — A bowl (Fig. 18) has the shape of an unknown sur- 
face of revolution; it seems probable from its form that it may have 
been generated by a parabolic curve. Find the time necessary for 
the level to fall from the height A2 to a height hi when three cross 




Fig. 18. — A bowl shaped like a surface of revolution generated by 
part of a parabolic curve. 

sections Fq, Fi, and F2 are given. Solution, — We assume that 
the siuf ace is generated by a parabohc curve of some unknown 
degree, and that the vertex of the curve lies at a distance z from 
the bottom of the bowl. The areas of horizontal cross sections vary 
then as some unknown power of the distance (z + x) from the 
vertex, and we have 

z + x= (aF)'' (7) 

This equation contains three unknown parameters, or determining 
constants of the surface, viz., 2, a, and n. These three may be de- 
termined by appljdng Eq. (7) to the three given cross sections. We 
have 

z=(aFor] 

z + hi=iaFir\ (8) 

z + h={aF2r] 
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Having determined the parameters of the surface by solving the 
three simultaneous Eqs. (8), we get from Eq. (7) 

p^(z_±x)}^ ^Qj 



Substituting this value of F in Eq. (3) and integrating, the re- 
quired time t is determined. 

In actual application to numerical problems this method presents 
two difficulties: (a) Eqs. (8) can be solved by trial only; (6) Eq. 
(3), with the expression (9) for F in it, can be integrated in finite 
form only if n is an integer; such a case would be a rare chance. 
Even then the expression is rather compUcated. If n is a fractional 
number, the integration has to be accomplished by Simpson's rule, 
or some such method. The procedure is illustrated in the next 
problem. 

Prob. 12. — Apply the method explained in the preceding prob- 
lem to the following data: Fq = 100 sq. cm.; Fi = 155 sq. cm.; 
F2 = 205 sq. cm.; hi = 10 cm.; hi = 26 cm. The area of the 
orifice is 2.5 sq. cm. Arts, t = 17.6 sec. 

Solution. — Subtracting the first of Eqs. (8) from the second gives 

10= (100a)«(1.55«-l); ..... (10) 
subtracting the first equation from the third, we get 
26= (100 a)** (2.05^-1). 

Eliminating (100 a)** from these equations by division gives an 
equation for n : 

2.6 X 1.55« - 2.05" = 1.6. 

To solve this equation, plot a curve giving the value of the right- 
hand side of the equation for different values of n; the required 
value of n is that for which the ordinate of the curve is equal to 1.6. 
Knowing n, the parameter a is computed from Eq. (10), and finally 
2 is determined from the first equation (8). 

In order to integrate Eq. (3) using Simpson's rule (see Chapter II), 
the height of the vessel between the levels and h^ is divided into 
an even number of equal parts, say 10 or 12, and for each hori- 
zontal plane of division the value of 

'-^. <"' 
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is calculated, using the value of F according to expression (9). 
Then, if AA is the distance between the planes of division, we have 



Jo Vx 3 



\yo + 4 (2/1 + 2/, + 2/, + 2/7 + 2/9 + 2/ii) 
+ 2 (2/2 + 2/4 + 2/6 + 2/8 + 2/10) +2/12} (12) 



After integration, Eq. (3) becomes 
cA V2g't= F, 



(13) 



from which expression the value of t can be calculated. The 
student is advised to perform these operations and to get the final 
numerical result. 

Prob. 13. — Show how to calculate the time necessary for empty- 
ing a reservoir of entirely irregular form, or the time necessary for 
lowering the level of the liquid by a giyen amount. This problem 
has an application to storage reservoirs, such as are formed by 
constructing a dam across a valley, for irrigation or power purposes. 
Hint, — Use equations similar to (11), (12), and (13) in the pre- 
ceding problem. 

Communicating Vessels. 

Two cylindrical vessels (Fig. 19) are connected by a 
short pipe of cross section A. The horizontal cross 
section of one vessel is Fi ; that of the other is F2. The 
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Fig. 19. — Two communicating vessels. 

original difference of the levels of the fluid is hi. It is 
required to calculate the length of time which it will 
take for the difference of levels to become h2. 
The orifice being submerged, the hydrauKc head is 
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equal to the instantaneous difference h of the levels in 
the vessels, so that, according to Eq. (3), 

cA V2gh 'dt= - Fidxi. . . . (14) 

The water which flows out of the first vessel raises the 
level in the second vessel, so that there exists the re- 
lationship 

-Fidxi=F2dx2 (15) 

The minus sign is necessary in this equation because Xi 
increases when X2 decreases. 

The variables Xi and X2 are also connected by the 
relation 

Xi — 0C2 = h (16) 

Eliminating Xi and Xa from these three equations, and 
integrating the result, the solution of the problem is 
obtained. We first differentiate Eq. (16) and elimi- 
nate dx2 by means of Eq. (15). The result is 

Substituting this expression in Eq. (14), we get 



cA V2gh.dt = -Fdh, . . . (17) 
where, for the sake of abbreviation, 

Eq. (18) is identical in form with Eq. (3), so that the 
result of integration is mathematically identical with 
Eq. (4), and we get 

cAV2g't = 2F{Vhi-Vh2). . . (19) 

This equation may be solved for t, which is the time 
required for reducing the difference of the levels from 
hi to h2. 
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Prob. 14. — Two rectangular tanks are given; the horizontal 
dimensions of one are 1.5 m. by 2.5 m., those of the other are 0.6 m. 
by 1 m. These tanks are connected by a pipe 5 cm. in diameter. 
The original difference of levels is 3 m. Plot a curve showing the 
time necessary to reach various differences of levels. 

Arts. h2 = 2 m.; t = 63 sees. 
A2 = 1 na.; t = 145 sees, 
fe = m.; t = 343 sees. 

Prob. 15. — What becomes of Eq. (18) when the area of one of the 
vessels is very large as compared to that of the other, as for instance 
in the case of a lock and a section of a canal considered above? 
Hint — If F2 is very large, divide the numerator and the denomi- 
nator by F2. 

Prob. 16. — Extend the theory of communicating vessels to the 
case when the vessels are not cylindrical, Fi being given as a func- 
tion of a:i, and F2 as a function of ^2. Apply the general method to 
the shapes of vessels given in Prob. 3, and in Probs. 7 to 12; also 
to various other shapes and their combinations. Hint, — Integrate 
Eq. (15) and eliminate X2 from the result by means of Eq. (16). 
Integrate Eq. (14) with respect to Xi instead of h, 

Infltix Simultaneous with Discharge. 

There are cases in practice when fluid is continu- 
ously added to a vessel, and at the same time some of 
it is being discharged. Such is the case in a storage 
reservoir, out of which water is being drawn for power 
or for irrigation; and at the same time there is an influx 
of water from a river above. Let q be the constant 
amount of Uquid added to reservoir I per unit of time 
(Fig. 20). Then Eq. (3) becomes: 

{cA y/2gx -q)dt= -Fdx. . . . (20) 

Since it makes no difference in what way the con- 
stant amount q is added, let us assume that it is being 
added from an infinitely large vessel II, through an 
identical orifice A, the level in reservoir II being main- 
tained constant and equal to m. This level m is deter- 
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mined by the condition that the discharge be equal 

to q; or 

cAV2gm = q (21) 

Equation (20) becomes then 

cA \/2g ( VS - Vm) dt= -Fdx. . . (22) 

Comparing Eq. (22) with Eq. (20) it will be seen that 
Eq. (22) is more sjmometrical, and therefore the solu- 




FiG. 20. — Influx simultaneous with discharge. 

tion can be more readily interpreted and analyzed. 
This is the reason for introducing the head m instead 
of the discharge q. 

If the horizontal area F of vessel I is constant, Eq. 
(22) may be integrated thus : 

cAV2^.t = F T—J^^. . . (23) 

In order to integrate the right-hand side of this equa- 
tion a new variable z is introduced: 

z = Vx — Vm, 

dx dx 



so that • 



dz = 



2Vx 2{z+ Vm) 
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Substituting, we get 
or, finally. 

Thus, Eq. (23) becomes 

cA V2i .t=2F\(Vh-Vh2) + Vmlog ^""^!!! T (24) 

Comparing this expression with Eq. (4) it will be seen 
that the last logarithmic term represents the effect of 
the influx in lengthening the necessary time of dis- 
charge. If there is no influx, m = 0, and Eqs. (24) and 
(4) become identical. 

Prob. 17. — Solve Prob. 1 when there is an influx of 5 liters per 
second. Can the tank be emptied, and if not what is the lowest 
level which the water can reach? 

Ans. A2 = 2.5 m.; q = 8.2 l./sec; t = sec. 
A2 = 1.25 m.; q = 5.84 l./sec; t = 798 sec. 
h2 = 0.92 m.; g = 5.0 l./sec; t = 00 sec. 
Prob. 18. — Investigate the value and the sign of t in Eq. (24) 
for hi> h2> m, h\> m> ^2, and m> hi> h^. 

Prob. 19. — Deduce an expression similar to (24) for the case of 
an inverted pyramid or cone with a small hole in the apex. Evalu- 
ate the integral f or w = 9 cm., from hi = 25 cm. to ^2 = 15 cm. 
Hint. — See Prob. 8 above. Ans, 2758 (cm.)^-^ 

Large Orifices. 

In the case of a large non-submerged orifice in the 
side of a vessel, the discharge of a fluid takes place 
under a different head at different parts of the orifice, 
so that Eq. (2) is applicable only to an infinitesimal hori- 
zontal strip of the orifice for which the head may be 
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considered constant. Thus, in the case of a rectangu- 
lar orifice of width Z, we have 



dQ = cdA 'V = C'ldh V2gh, . . . (25) 
and the total discharge 

' Vhdh, 

where A^i and fe are the distances from the free surface 
of the liquid to the upper and the lower edge of the 
orifice respectively. Integrating, we obtain 

Q = f cZ (fe V2^ ~ ki V2gki), . . (26) 

also sometimes written in the form 

Q = fcZV^(fc2*-M). . . . (27) 

Form (26) is perhaps more convenient for computa- 
tions, because the values of V2gk may be taken directly 
from available tables. 

Prob. 20. — A rectangular tank 15 m. by 60 m. is filled with 
water to a level of 3.2 m. above the bottom. How long will it take 
io empty half of the tank through a rectangular orifice 0.5 m. wide 
and 0.2 m. high, the lower edge of the orifice being 0.1 m. above the 
bottom? 

Solvtion, — By analogy with Eqs. (2) and (3), and using expres- 
sion (26), we get for an infinitesimal element of time dt corresponding 
to a change dx in the level: ' 



J X 0.6 X 50 [{x + 20) V2g(x + 20) - x V2gx]dt 
= - 1500 XGOOOdXy 



or, separating the variables, and integrating, 

2 2tXlO~^= —-= l^ 

V2fir Ji3o {x + 20) V(a; + 20) -xVi' 

The integral can be computed readily by using Simpson's rule. 

Ana, t = 376 sees. 
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Prob. 21. — Deduce a formula corresponding to Eq. (26) for a 
triangular orifice (Fig. 21), the apex being downward. Hint, — 
In Eq. (25) Z is a linear function of h. 

Am. Q = Octane V2g[2k2^-'5kjci^ + dki^l 




h-i- 



FiG. 21. — A large triangular orifice. 

Prob. 22. — Deduce a formula corresponding to Eq. (26) for a 

parabolic orifice. See Fig. 22, in which y = px^, 

Am. Q = \Tcy/2gk^/ V^, 



|<— <x_^ 





Fig. 22. — A parabolic orifice. 



Fig. 23. — A round orifice. 



Prob. 23. — Deduce a formula corresponding to Eq. (26) for a 
round orifice. Solviion. — With the notation according to Fig. 23, 

we have 

1 = 2 V;^^2^ 

so that Eq. (25) becomes 

dQ = 2c Vr2 - 2/2 . V2g{h -f y) •dy. . . . (28) 

This expression must be integrated from ?/= —rtoy= +r. Since 
it cannot be integrated in a finite form, two ways are open: 
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(1) If only a numerical result is desired, use Simpson's rule. 

(2) If a general expression for Q is required, expand expression 
(28) into a convergent infinite series and integrate each term sepa- 
rately. Of course, the result is of use only if it also happens to be 
a convergent series. 

In order to make expression (28) integrable we expand ^h + y 
according to the binomial theorem, so that 

dQ = 2cV2^.Vj^3^jl + ||-|g+ . . .j%. 

We thus have, using a table of integrals: 

j^y/j^-y^*dy = i^ Vr^ - y^ + r^sin-i (2//r)J_^= Jirr*; 

j^^y V?~:^^dy = - i[V(r2-t/2)']^= 0; 

£^V V'y^^^^y = [-| Vi^^2+|'(2/v'i^^^+r^sin-i^^]^^ =^, 



The first of the foregoing integrals represents the area of the half 
circle (left- or right-hand side of the orifice); the second integral 
represents its static moment, and the third integral its moment of 
inertia with respect to the horizontal axis passing through the 
center. Thus, in this case the results could have been written 
without actually performing the integration. It will also be seen 
that all the terms containing odd powers of y vanish in the 
result. 

Thus, finally, 

or 

= c^^.,r^jl-l(^J j. ... (29) 

Without the correction terms in the parentheses, the formula is 
identical with Eq. (2) for a small orifice. 

Prob. 24. — Compute two more terms in Eq. (29). 
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Prob. 25. — The discharge through various large round orifices 
is to be computed using the approximate formula (2), Plot a curve 
giving per cent error in the result for different values of r/h, 

Prob. 26. — The diameter of a round orifice is 1.4 m., the dis- 
tance from the upper edge to the level of the water is 0.3 m. Com- 
pute the discharge, using Simpson's rule, and compare the result 
with that obtained from Eq. (29)* 
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CHAPTER IV. 

FORM OF LIQUID JET. 

A JET of liquid AC (Fig. 24) is issuing through a 
horizontal nozzle A in the side of the vessel MN. 
The level of the liquid in the vessel is kept constant. 
The problem is to find the mathematical shape of the 
jet, knowing its initial velocity v. 




Fig. 24. — A jet issuing from a horizontal nozzle. 

The motion of a particle of liquid, such as P, may be 
resolved into a uniform horizontal motion at the initial 
velocity due to inertia and into an accelerated motion 
downward imder the action of gravity. If x and y 
are the coordinates of the point P, referred to point A 
as the origin, the foregoing components of the motion 
are expressed mathematically thus: 

x = vt, (1) 

43 



(2) 
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where t is the time, counted from the moment at which 
the particle P issued from the nozzle,, and g is the 
acceleration due to gravity. 

EUminating t from Eqs. (1) and (2) the required equa- 
tion of the path is obtained. Before doing this, Eq. (2) 
must be integrated. The first integral is 



dt 



= gt + Ci, 



where Ci is the constant of integration, which is deter- 
mined by the specific conditions of the problem. 
Naniely, when t = 0, the vertical component dy/dt of 
the velocity is also zero, since the nozzle is horizontal. 
Hence Ci = and the first integral is 

Integrating again, we get 

When t = 0, the ordinate i/ = 0, so that C2 = 0, and 
we finally have as the integral of Eq. (2) 

y = h9t' (3) 

Eq. (3) is sometimes deduced in an elementary way 
as follows: The acceleration g being constant, and the 
initial velocity zero, the final velocity after an interval 
of time t is gt, because, by definition, acceleration is 
the rate of increase of velocity per unit of time. Thus, 
the average velocity during the period of time t is ^ gt, 
and hence the length of the path in the direction of g 
IS y = 2 9^ ' i == ^ 9^' Relation (3) may also be con- 
sidered as a direct result of experiments on bodies fall- 
ing under the action of gravity. 
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Eliminating t from Eqs. (1) and (3) we obtain 

y-i;^'^ (4) 

This is a parabola with its vertex in the origin. 

Prob. I. — Show that Eq. (4) is physically homogeneous. In 
other words show that the dimension of the right-hand member is; 
that of a length which is also the dimension of y. Solviion, — 

(length/tmae)2 ^ ° ' 

Prob. 2. — A jet issues from the vessel MN (Fig. 24) at a velocity 
of 3 m. per sec. and is delivered into the vessel c, 1.25 m. below the 
nozzle. Plot the shape of the jet and determine the horizontal 
distance between the nozzle and vessel c. Arts, x = 1.515 m. 

Prob. 3. — Show that when a discharge takes place imder ideal 
conditions, that is, without friction or eddies in the nozzle, the 
directrix of the parabola lies in the free surface QN of the liquid. 
Solution, — When there are no frictional resistances to the flow 
of Uquid, the total potential energy of a particle of mass m on the 
free surface is converted into kinetic energy as it issues from the 
nozzle. Thus, if A is the head, we may write: (ffm) h = i irn^f 

hence 

1;= V2gh (5) 

Expression (5) is called Torricelli's theorem, and it may be con- 
sidered as a well-established experimental fact. Substituting in 
Eq. (4) the value of v from Eq. (5), we obtain 

x2 = 4% (6) 

Comparing this with the standard form of the equation of a pa- 
rabola, x^ = 2pyy and remembering that the distance between the 
vertex and the directrix is 1 p, the required relation is proved. 

Prob. 4. — Prove that when Eq. (6) holds true, the range DC 
of the jet (Fig. 24) is equal to twice the abscissa AK of a semi- 
circle described on ND as a diameter. 

Prob. 5. — Prove on the basis of the result of Prob. 4 that to 
obtain a maximmn range the nozzle A must be placed at such a 
distance h from the free surface as to divide ND in half; also that 
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the range is the same for any two positions of A equidistant from 
N and 2). 

Because of frictional resistances and of eddies, the 
actual velocity of discharge is lower than that deter- 
mined from Eq. (5). Therefore, in practice, this 

equation is modified to 

v=fV2^, (7) 

where / is the so-called coefficient of velocity.* If, 
for instance, / = 0.95, the actual velocity is equal to 
95 per cent of the theoretical velocity corresponding to 
the same head h. Five per cent of velocity is lost in 
friction and eddies in the nozzle. The value of / de- 
pends upon the shape and size of the orifice and of the 
nozzle, and on some other conditions of discharge. 
Introducing in Eq. (4) the value of v from Eq. (7), 
we get 

x2 = 4%.A (8) 

BO that the range x decreases with the coefficient /, and 
one is a direct measure of the other. This relationship 
provides a simple means for an experimental determi- 
nation of /. 

Prob. 6. — A horizontal jet (Fig. 25) issues from the nozzle A, 
the center of which is hi m. below the free surface; another jet 
issues from the nozzle B which is ^ m. directly below the first. 
Find the coordinates of the point at which the two jets meet each 
other, assuming the coefficient of velocity of the two nozzles to be 
the same. Ans. Xi = ^4:hifi^{hi + h). 

Prob. 7. — Solve Prob. 6 when the coefficients of velocity are 
different for the two jets. Can the coefficient of velocity of the 
lower jet be made so small that the jets will not intersect? 



Am. ^^-\f^^kr + h)-kf^' 



* This coefficient is different from the coefficient of discharge c, used 
in the preceding chapter. 
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Prob. 8. — The coordinates of the following points have been 
actually measured on a very thin jet (Fig. 24) : 

y = I 2 3 meters. 

X = 2.77 3.92 4.80 meters.- 

The head h equals 3 m. Determine the coeflSicient / of the nozzle. 
Hint. — Mark the given points on a sheet of logarithmic cross- 
section paper and draw the nearest straight line through them at 
a slope of 2 : 1. Am. f = 0.8. 

The foregoing theory may be generaUzed for the case 
when the axis of the nozzle is inclined at an angle a 




Fig. 25. — Two intersecting jets. 

to the horizontal (Fig. 26). The horizontal projection 
of the initial velocity v is now v cos a, so that Eq. (1) 
becomes 

X = vcosa^t (9) 

The vertical projection of the initial velocity is —t; sin a, 
so that Eq. (3) becomes 

y^^gfi — V sin a»t (10) 

This is because dming the time t a particle not only 
falls imder the action of gravity by the amoimt ^ gf, 
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but it is also carried by inertia upward on accoimt of 
the vertical component t; sin a of the initial velocity t;. 
EUminating t from Eqs. (9) and (10) the following equa- 
tion of the shape of the jet is obtained: 



y = 



2 V^ COS^ a 



xtBJla, 



(11) 



For a horizontal nozzle a = and Eq. (11) becomes 
identical with Eq. (4). 




Fig. 26. — Jet with an inclined axis. 

Introducing the value of v from Eq. (7), as before, 
Eq. (11) becomes 



y = 



4 hf cos^ a 



— X tan a. 



(12) 



Prob. 9.* — Eq. (12) represents a parabola; reduce its equation 
to the standard form x^ = 2py, Solution, — Let xq and yo (Fig. 26) 
be the unknown coordinates of the vertex K of the parabola. The 
equation of the parabola is of the form 

(x-Xoy = 2p(y-yQ), 
or x^ + xo* — 2 xa^o = 2 py — 2 pyo. 

* Most of the problems that follow apply also to the finding of trajec- 
tories and ranges of projectiles, except that the resistance of the air and 
other disturbing influences modify somewhat the path and the velocity 
of a projectile. 
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The parabola passes through the origin A, so that this equation 
must be satisfied when a; = and 2/ = 0. Substituting these values 
we get the relation 

Xo^=-'2pyo, (13) 

so that the foregoing equation is simplified to 

x^-2xxo = 2py (14) 

Comparing it term by term with Eq. (12) we find that 

p = 2hPG08^a, (15) 

Xo = 2hP cos2 a . tan a, .... (16) 
and using Eq. (13) 

yo= -hfcos^ ataxia (17) 

The coordinates xq and yo are expressed through the parameter 
2p thus: 

^« = ^^^: , i (18) 

2/0 = -iptan^a) 

Prob. 10. — A jet issues at an angle of 30 degrees to the hori- 
zontal and at a velocity of 15 m. per sec. Find the origin of the 
parabola and construct the curve using the value of its parameter 
2 p, Ans. xo = 9.93 m.; 2/0 = —2.87 m. 

Prob. II. — Find the value of xq directly from Eq. (12), as the 
point at which the height of the jet reaches its maximum. Check 
the value so found with that given by Eq. (16). 

Prob. 12. — How would you find from Eq. (12) the horizontal 
distance at which a vessel C must be placed in order to receive the 
liquid? Explain the meaning of the double sign in the solution. 

Prob. 13. — Find the position of the directrix of the parabola 
in Fig. 26. Solution. — The height of the directrix above the ver- 
tex is J p, so that from Eqs. (15) and (17) its height above the 
nozzle is hp cos^ a tan* a + hp cos* a = hp. Hence the directrix al- 
ways lies below the free surface of the liquid, and coincides with it 
in the ideal case of / = 1.00. 

Prob. 14. — Find the value of angle a corresponding to the 
greatest possible horizontal range AL at the level of the orifice. 
Solution, — In Eq. (12) let y = 0; then, factoring out x/cos a, the 
result is x/(4 hp cos a) = sin a, or a; = 2 hp sin 2 a. The abscissa 
a; is a maximum when sin 2 a = 1 ; thus, Xm&x = 2 hp, when a = 45°. 

Prob. 15. — Water flows out of a vessel through an inclined 
nozzle, at a head h equal to 3 m. Two ordinates of the jet have been 
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measured: 2/1 = at a:i = 3 m. and y^ = 3.5 m. at 0:2 = 6 m. 
From these data calculate angle a of the inclination of the nozzle, 
the coeflSicient of velocity /, and the highest point K which the jet 
reaches above the nozzle. Am. f = 0.758; a = 30° 15'. 

Prob. 16. — Solve Probs. 6 and 7 when both nozzles are inclined 
at the same angle a. / — , . ., — —rr — 

"■"■ '-»«— v^i^^- 

Prob. 17. — Solve Probs. 6 and 7 when the angles ai and az of 
inclination of the jets are different. 

Prob. 18. — The following three points were obtained experi- 
mentally on a jet at A = 3 m. : 

x= 2 4 6 m. 

y = -0.46 + 0.468 + 2.884 m. 

Construct by trial a parabola which comes the nearest to these 
points, and from it determine the angle a of inclination of the nozzle 
and the velocity coefficient /. Note. — For a more accurate method 
based upon the method of least squares see Prob. 31 below. 

Am. a = 30°; / = 0.80. 
Prob. 19. — Find angle a at which a jet must issue from a fire 
hose at a given velocity in order to strike a given point. Solution. — 
Eq. (11) may be reduced to the form 

a sin 2 a + 6 cos 2 a = A;, (19) 

where a, 6, and k are known quantities. This equation is solved 
by introducing an auxiliary angle <t> such that 

6/0 = tan (20) 

Substituting, Eq. (19) becomes sin 2 a + tan cos 2 a ^ k/a, or 
sin (2 « + 4) _ k 
C08<t> a' 

From this equation the unknown angle (2 « + <t>) may be cal- 
culated and then a is determined, knowing <t> from Eq. (20). 

Prob. 20. — A jet issues from a vessel under the head ^ = 3 m. 
Find the angle of inclination of the nozzle, knowing that the jet 
strikes a point whose coordinates are x = 2.5 m. and y = 1.35 m. 
Coefficient of velocity/ = 0.90. 

Am. a = 76° 27', or a = - 14° 50'. 

Prob. 21. — A jet whose shape is represented by Eq. (12) strikes 
an inclined floor. Find angle a of the nozzle at which the range on 
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the floor is a maximum. The position of the floor is given by the 
following data: Its plane is perpendicular to that of the jet; it is 
b m. below the nozzle, and it rises at a slope of 1 ; n. 

First SolvUon, — The equation of the trace of the floor in the 
plane of the jet is 

2/ = ^-^ (21) 

n 

Combining this with Eq. (12) in order to find the point of inter- 
section of the two curves, we get 

^--=n7r^--«tana (22) 

To find a maximum range this equation must be solved for x, 
and then a determined from the condition 

^=0 (23) 

Oct 

Second Solution. — The preceding solution leads to rather long 
formulaB; therefore, it is advisable to introduce condition (23) 
before solving Eq. (22) for x. Dififerentiate Eq. (22) with respect 
to a, considering x as a function of a: 

_ 1 (^\ = 2 X (dx/da) . 2 x^ sin a _ . /dx\ _ x 

n\da) AhPcOS^a 4^/2 cos' a \da/ COS^ a 

In view of condition (23), this equation is simplified to 
2 x^ sin « _. X 

4^/2 cos' a COS^a' 

or 

xtajia = 2hp (24) 

Solving this equation with Eq. (22) as two simultaneous equations, 
the required values of x and « are obtained. The ordinate y is then 
determined from Eq. (21). 

Ans, 2 + b/(hP) = 2/(n tan a) + l/sin^ «. 
Prob. 22. — A jet A (Fig. 27) is used for watering a garden sit- 
uated on a slope inclined to the horizon at an angle /8. Find the 
angle a at which the jet must be held in order to reach the furthest 
possible point B, The jet is supposed to issue directly from the 
ground, at a given velocity v. Solution. — The equation of the 

straight line AB is 

y=—x tan/8 (25) 
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Solving it with Eq. (11) in order to find the point of intersection B, 
we find 

sm (a — /8) . 



X = 



The maximum range is found from the condition 
dx 2v^ 



I cos (a— j8) cos a — sin a sin (a — /8) I = 0, 



da ^ cos /8 I 
or cos (2 a - /3) = 0. Hence, 2 a - /3 = 90°, or 

a= i/3 + 45°. . . . 



(26) 




Fig. 27. — A jet which strikes an incUned plane. 

Prob. 23. — At what angle a and at what velocity v must the 
jet be issued (Fig. 27) in order to strike a given point 5 at a given 
angle <t> to the horizon? Solution. — Let a be the horizontal dis» 
tance between A and B, We have from Eq. (11) 



or 



■f— tan a = tan ^. 



(27) 



Applying Eq. (11) to point B gives 
a tan j8 = 



2t;2cos2a 



— a tan a, 



or 



-T^ 2 tan a = 2 tan /8. 

t^ COS^ a 



(28) 



Subtracting Eq. (28) from Eq. (27), tan a is determined thus: 
tan a = tan — 2 tan /8. 



Digitized by 



Google 



Chap. IVJ FORM OF LIQUID JET. 53 

Prob. 24. — A fountain jet issues through a nozzle (Fig. 28) 
provided with small openings at the end, which throw water in all 
directions. Find the equation of the surface enoeUrping all the 
individual jets. 

First Solution, — Since the required surface is a surface of revolu- 
tion with a vertical axis, it is sufficient to consider the curve of its 




Fig. 28. — A vertical nozzle with small holes at the end. 

intersection with the plane of the paper. Individual jets are repre- 
sented by Eq. (12), angle a being different for each jet. It will be 
seen from the figure that the enveloping curve touches each jet 
at the point of its maximum range (Prob. 23). Hence, relation 
(24) holds true for each point common to a jet and to the envelope. 
Therefore, eliminating a from Eqs. (12) and (24), we obtain a gen- 
eral relationship between the coordinates of all the points of tan- 
gency, in other words, the equation of the envelope. We have 

from Eq. (24) 

1 4 A^f* 

^ =l + tan2cr=l + ^^-^ 



COS^a X^ 

and substituting this value into Eq. (12), we get, after reduction, 
x^ = 4hP(y + hP), (29) 
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which is a parabola. Hence, the required surface is a paraboloid 
of revolution. 

Second Solvtion, — ^ According to the general rule for finding an 
envelope, let 

/(^,2/,«) = (30) 

be the general expression for a family of curves (in our case, all the 
jets in one plane), where a is the variable parameter which dis- 
tinguishes each curve; in our case it is an angle. The equation of 
the envelope is obtained by eliminating a from Eq. (30) and from 
the derivative equation 

^My;«) ^0 (31) 

da 

In our case Eq. (30) corresponds to Eq. (12), and its partial deriv- 
ative with respect to a is Eq. (24). Thus, we have the same solu- 
tion as before. 

Prob. 25. — Show that the enveloping parabola, Eq. (29), has 
its vertex in the plane of the " efifective " level of the liquid in the 
supply tank, and that the focus of the parabola coincides with the 
center of the hemispherical end of the nozzle. Note, — The efifec- 
tive level corresponds to the head hp which would be maintained 
for the same velocity of the jet, if the discharge were ideal, with- 
out friction. 

Prob. 26. — Using the data of Prob. 10, plot a curve the ordinates 
of which will give values of velocity at different points of the jet. 
Also plot another curve the ordinates of which will give values of 
the trigonometric tangent of the angle which the jet makes with the 
horizontal at those points. Plot these curves to the same abscissae 
as the jet itself. Hint. — Velocity = V{dx/dty + (dy/dty; use 
Eqs. (9) and ( 10), and eliminate L 

Am. V = V225- 11.33^ + 0.57x2; dy/dx = (0.755 x - 7.5)/13. 

Prob. 27. — Show that in the preceding problem the curve of 
tangents is a straight line, and the curve of velocities is a hyperbola 
the vertex of which is on the abscissa at which the jet reaches its 
highest point. 

Prob. 28. — Find the length of the jet, along the curve, expressed 
by Eq. (12), from the nozzle to the highest point. 

Am. s=hp\sm a+cos^ a loge (1 +sin a) — cos^ a log. (1 — sina) } . 

Prob. 29. — Find the angle of inclination of the nozzle such that 
the area included between the jet and the horizontal line through 
the nozzle will be a maximum. 
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First Solvtion, — Referring to Fig. 26, area AKL must be a 
maximum. For a parabola this area is equal to } of the area of the 
circumscribed rectangle, so that the product Xol/o must be a max- 
imum. Therefore, according to Eqs. (16) and (17), the product 
COS^ a tan' a must be a mflYiTn unn^ or 

cos a sin' a = max (32) 

Taking a derivative and equating it to zero we find that sin « = 
i V3, or a = 60°. 

Second Solution. — Integrating Eq. (12) directly, we find that 
the area AKM is 

^0 ^ 12hPc08^a 2 

Substituting the value of xq from Eq. (16) gives A = Const. X cos a 
sin' a, so that the result is identical with that obtained from 
Eq. (32). 

Prob. 30. — Solve Prob. 8 using the method of least squares. 
Solution. — Let the unknown (most probable) parameter of the 
parabola be 2 p, and let some approximate value of it be 2 P. This 
value may be obtained by a rough trial and substitution into one 
or two sets of the experimental data. Let p be the unknown cor- 
rection, so that 

2p = 2P + p, (33) 

and the equation of the most probable parabola becomes 

x'2=(2P + p)2/' (34) 

The sign " prime " means that x and y in this equation are not the 
actually observed values, but the most probable values. Since the 
observed points presimiably do not lie exactly on this parabola, 
we have 

x^ — {2P + p)y = error, 

where x and y refer now to the actual experimental data. For the 
most probable parabola the sum of the squares of these errors must 
be a minimum, or 

S[a:2-(2P + p)2/P = min. 

While a long detailed discussion of this proposition would be out 
of place here, the student can readily see that the condition 

S errors = min. 
would not give the most probable curve. The individual errors 
may be very large, but some positive errors may cancel negative 
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errors making the sum quite small. On the other hand, the con- 
dition 

S (error)2 = min. 

necessitates thai the individual errors be small because all the 
terms in the summation are positive.* 

To find a minimum of the preceding expression, we equate to 
zero its derivative with respect to p, and obtain 

S2/[x2-(2P + p)2/] = 0, 
or 

from which 

^, sy(yp,) (3,) 

Since 2 P is known, small corrections 

rj' = x^-2Py (36) 

for the experimentally observed points can be calculated, and then 
p computed from the formula 

''X^ (^^) 

The student is advised actually to perform these computations, 
and to check the answer to Prob. 8, because only in this way can he 
fully appreciate the advantages of the method of least squares. 

Prob. 31. — Solve Prob. 18 using the method of least squares. 
Solution. — According to Eq. (14), the most probable parabola is 
determined by the parameter 2 p and the abscissa xq of the vertex. 
Let 2 P and Xo be approximate values of these quantities deter- 
mined by trial, and let p and ^ be the unknown corrections. Then 
the most probable values are 

2p = 2P + p) 
Xo = Xo + ^ 
and Eq. (14) becomes 

x'^-2x'(X + ^) = (2P + p)y\ 



'[ (38) 



where a;' and y' are the coordinates of the most probable parabola. 
If, however, instead of these, we substitute the coordinates x and y 

* See also the last portion of Chapter III, in the part on "Thermo- 
dynamics. " 
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from the observed curve, the equation ceases to be satisfied, and we 
have 

error = x' - 2x (X + ^) - (2P + p)y. 

For the most probable values of ^ and p the condition must be satis- 
fied that 

S (error)* = S [x* - 2a; (X + ^) - (2P + p) yp =niin. 

Taking partial derivatives with respect to ^ and p, and equating 
them to zero, we obtain the following two equations: 

Sx [x2 - 2 a: (Z + ^) - (2 P + p) 2/] = ; 

S2/[x2-2a:(X + ^)-(2P + p)2/] = 0; 
or 

2 ^2x2 + pXxy = Xx {x^ -2xX- 2Py); 

2^i:xy + pS2/2 = Xy (x* - 2xX -2Py), 

In these expressions we denote the discrepancies by ly*, as in the 
preceding problem, so that for any experimental point rf =^ x^ — 
2 xX — 2 Py, where 17' is a small quantity. Solving the two pre- 
ceding equations as simultaneous equations for p and ^, we obtain 



^ ' Sx^Sy* - (2x2/)2 



Sa;2S2/2 - {l^xyY 



(39) 



Having determined p and ^, the quantities 2 p and xo become 
known from Eqs. (38). Finally, angle a and coefficient / may be 
calculated from Eqs. (15) and (16). The student is advised to per- 
form these computations in order to become familiar with the 
method of least squares. He will probably find that the accuracy 
of an ordinary slide-rule is not sufficient for such computations, and 
that either a more accurate computing device or a table of log- 
arithms must be used. 
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CHAPTER V. 

PROBLEMS ON WATER SUPPLY. 

The usual method of supplying water to cities con- 
sists in pumping it from an available source, such as 
a river or springs, to an elevated reservoir. From the 
reservoir the water is fed through pipes by gravity 
to various parts of the city. The level of the water in 
the reservoir must be high enough to supply the neces- 
sary pressure at the highest points in the distributing 
pipes, and also to overcome a considerable loss of 
pressure in the pipes themselves, due to frictional re- 
sistance. Therefore, in the solution of problems in- 
volving diameters of pipes it is necessary to know in 
what way the loss of pressure depends upon the diam- 
eter of a pipe and upon the average velocity of water. 
The fundamental formula for the loss of pressure per 



hi 



H 



! ! 

Fig. 29. — Fall of pressure along a pipe. 

imit length of a pipe will be f oimd in almost any stand- 
ard text-book on hydrauUcs, but for the convenience 
of the student it is deduced below. 

Let AB (Fig. 29) be a portion of a pipe of length I 
and diameter D, and let water flow through this pipe 

58 
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at a rate of Q cubic units per second. The free pressure 
of water at different points along the pipe may be ascer- 
tained by inserting small vertical tubes, called piezom- 
eters, shown at A and B. When the water in the pipe 
is at rest, the level of water in these small pressure 
coliunns rises to that in the supply tank. When the 
water in the main pipe is in motion these levels sink,, 
on account of loss of pressure due to frictional resistance 
in the pipe between the tank and the section of the 
pipe under consideration. Thus the difference of the 
levels at A and at B is due to the frictional resistance 
along AB. 

The pressure of the water on cross section No. 1 of 
the pipe is ^ tt D^hiW, where w is the weight of a cubic 
unit of water. Similarly, the pressure on cross section 
No. 2 is ^ TT D^JhWy so that the net force moving the 
volume of Uquid between A and B is ^ tt Dhv {hi - h^). 

Since the motion of the water is uniform, this force 
is balanced by the force of friction between the pipe 
and the water. This friction may be expressed as 
f 'wDly where / is the force of friction per unit area, 
and TT Dl is the inner surface of the pipe. Thus we have 

iTDhoihi-h) =f'TrDl. . . . (1) 

Friction / is found by experiment to be approximately 
proportional to the square of the velocity of the hquid, 
or 

f^av', (2) 

where a is an experimental coefficient. But the aver- 
age velocity v and the discharge Q through the pipe are 
related as follows: 
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SO that 

/ _ 16 gQ^ .o^ 

^^'^^^ ^^^ 

Substituting this value of / in Eq. (1) and combining 
all the constants, we obtain the following formula for 
the loss in pressure in terms of the difference in the 
levels of the piezometer tubes: 

hi-h2 = dQyD\ 

where c is an empirical coefficient. 

In practice, it is often convenient to know the loss 
in pressure per unit length of a pipe. Dividing both 
sides of the preceding equation by I we get 

i = {hi - h)/l = cQ^/D' (4) 

In this formula i is known as the hydraulic gradient. 
In Fig. 29, i is equal to the trigonometric tangent of the 
angle which KL makes with the horizontal. 

The coefficient c is not constant even for a given 
quaUty or roughness of the inner surface of the pipe, 
but depends somewhat upon the diameter of the pipe 
and also upon the velocity of the water. For our pur- 
poses c will be assiuned constant and equal to 23.6, 
when Q is expressed in Uters per second and D in centi- 
meters. The gradient i is a numeric, but both the loss 
of head (hi — h^) and the length of pipe I must be ex- 
pressed in the same units, for instance in meters. 

There are various tables in existence which give more 
accurate values for i than does formula (4) ; however, 
in order to treat certain problems mathematically, the 
value of i must be expressed through a formula. Then, 
after the solution has been found, it may be corrected 
by using a more accurate value of the constant c, taken 
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from a table, to correspond to the required diameter 
of the pipe and to the actual velocity of flow (see also 
Prob. 4 below). 

Prob. I. — What is the loss of pressure in a 10 cm. pipe 2.5 
kilometers long, for rates of discharge varying from zero to 30 liters 
per second; plot a curve. Ans, hi — hi = 59 Q^ cm. 

Prob. 2. — It is desired to represent the loss of head, according 
to formula (4), in the form of a chart of curves from which the loss 
of head per 100 m. of length of pipe could be read ofif directly 
for pipes ranging from 5 to 15 cm. in diameter in steps of one centi- 
meter, and for rates of discharge such that the velocity of flow would 
be between 0.5 and 2 m. per sec. What is the best way to arrange 
such a chart? Calculate a few points, and show the general char- 
acter and the limits of the curves. 

Prob. 3. — The relationship between i, Q, and D in Eq. (4) corre- 
sponds to a surface, if i, Q, and D be plotted along three axes of 
coordinates in space. Describe as accurately as possible, both in 
mathematical and in physical terms, the general character of the sur- 
face. Illustrate the description by free-hand sketches, showing a few 
cross sections of the surface. 

Prob. 4. — A formula for the hydraulic gradient, more accurate 
than formula (4), is 

(5) 






where A; is a numerical coefficient, and m and n are exponents some- 
what different from 2 and 5. Determine the values of k, rriy and n 
from the results of the experiments given in Fig. 30. Hint, — 

log z = log A; + m log Q — n log D (6) 

The experimental curves being replotted to a logarithmic scale give 
a series of approximately parallel straight lines. The tangent of 
the angle of inclination of these lines with respect to the Q axis gives 
m, while n may be computed from the horizontal distance between 
the lines. Arts, k = 22; m = 2.5; n = 4.8. 

Prob. 5. — Solve Prob. 4 using the method of least squares. 
Solution. — The method is explained in Probs. 30 and 31 at the end 
of Chapter IV. For the sake of uniformity it is convenient to repre- 
sent formula (5) in the form i = 10*'Q"*/D'*, because it can then be 
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reduced to r + m log Q — n log Z) = log i, with r, m, and n as un- 
known quantities. Let Af, AT, and i? be some approximate values 
of m, n, and r, determined by trial, and let the most probable or 
the best values be 

W = ilf + /i, 

n=^ N +Vj 
r =R +p, 
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Fig. 30. — Experimental curves of loss of pressure in certain pipes. 

where /*, j/, and p are unknown small corrections. Hence, we ob- 
tain p + /Li log Q — J/ log Z) = 17, where 17 is a known quantity. The 
condition for the most probable values is 

2(p + /ilogQ - J/logZ) - ,7)2 = min., 

and the problem is reduced to the solution of three linear simulta- 
neous equations. iVofe. — A seemingly shorter way is not to find 
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any approximate values, M, N, and i2, but to apply the condition 
of least squares to the exponents r, m, and n directly. In this case 
error = r + m log Q — n log D — log i, 



and 



2 (error)2 = S (r + m log Q — n log Z) — log iy = min. 



The three partial derivatives of this expression, equated to zero, 
furnish three simultaneous equations for the unknown exponents. 
A serious disadvantage of such a procedure is that it leads to very 
large numbers with several significant places, which must be care- 
fully computed in order to obtain a result that is at all accurate. 
This is because the result depends upon differences of large expres- 
sions almost equal to each other, so that the last significant figures 
really determine the result. On the contrary, with the method out- 
lined before, the simultaneous equations are appUed only to small 
corrections of the quantities already known, so that a considerable 
error in the computations affects the final result but Httle. A 
thoughtful student will derive considerable benefit from solving the 
problem both ways. 

Formula (4) applies to the simplest case when the 
discharge Q is constant along the pipe. It is not appli- 
cable directly to a street main, from which water is 
taken by house pipes, because the discharge is variable 



Hf 



P4R. 



dic^ I*- 



=pl+R 



'aip+R 



Fig. 31. — Loss of head in a pipe with a continuous uniform loss 
of discharge along its length. 

along such a main. While water is taken at a finite 
number of points, it is mathematically simpler to as- 
sume the loss of discharge to be continuous along the 
main. In this case formula (4) may be apphed to an 
infinitesimal length dx of the pipe (Fig. 31), and the 
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result integrated over the desired length of the pipe. 
Branch pipes to individual houses occur in practice 
at somewhat irregular intervals, and the demand for 
water is never the same or simultaneous from house to 
house. However, to simplify the conditions, the great- 
est demand for water is assmned to be continuous, and 
this demand is assumed to be uniformly distributed 
along the street main as if it were provided with a 
narrow longitudinal slot through which a uniform leak- 
age occurs. 

Let this side-demand be p Uters per second per Unear 
meter of the street. Let the pipe also carry an addi- 
tional quantity R of water, which is fed at the end of 
the pipe into the next street main. 

The quantity of water carried by the pipe at a dis- 
tance X from the end is 

Q, = R+xp (7) 

Substituting this value in formula (4) we get for the 
drop dH of pressure along the element dx 

dH = idx = c{R+ xvY dx/D\ . . (8) 

Integrating this expression over the total length of 
the main we find 

H = c f\R + xvY dx/D^ = c{lR^ + V^Rp + \ IY)/DK 

Let the total demand along AB be denoted by P, so 
that P = pi; the above expression reduces to 

H = cl{R'+RP + lP^)/D', ... (9) 

which gives the required relation between the loss of 
pressure H, the diameter D of the pipe, and the dis- 
charges R and P. 
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Prob. 6. — A 30 cm. main 2 km. long must supply 70 liters per 
second to the street as a whole. How much may it discharge into 
the next street main provided that the total drop in pressure should 
not exceed 2 kg. per sq. cm.? Ans. 64.3 l./sec. 

Prob. 7. — Show that when i^ = the loss of head is one third 
of what it would be if P were concentrated at B instead of being 
imiformly distributed along AB, 

Prob. 8. — Deduce the equation of the curve ab of hydraulic 
gradient (Fig. 31). Select the axes of coordinates so as to reduce 
the curve to its simplest form. Hint. — Integrate Eq. (8) between 
and x; the curve is a cubic parabola. 

. Ans. y=(jj^ofi-R^)/3p. 

Prob. 9. — Deduce an equation analogous to Eq. (9) when the 
distributed demand is on part k of the pipe only (Fig. 32). Check 



1 



P+R 



V 



< 12 * 



^ 



Fig. 32. — Loss of pressure in a pipe with a uniform demand over 
part h of its length. 

the answer for the cases previously considered when fe = 0, and 
when li = lz = 0. 

Ans. ^=;^5{ ^^^f^"^^'^ +^^(2?i + ^)+^^+^(^i + ^+MJ- 

Prob. 10. — Eq. (9) becomes identical with Eq. (4) if an equiva- 
lent demand Q is introduced, such that Q^ = R^ + PR + J jP. 
Show that 

R + ^>Q>R+iP, 
V3 

so that the usual assumption Q = R + 0.55 P is justifiable in 
practice. 
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Prob. II. — The character of a street is such that the demand p 
per unit length cannot be assumed constant. This demand is equal 
to po per unit length at A and is equal to pi at ^; p may be assmned 
to vary according to the straight-line law between these points. 
How is formula (9) modified? 

D\ ^3 3 ^ 20 ^ 20 15 / 

Prob. 12. — Modify and integrate Eq. (8) using for hydraulic 
gradient, expression (5) instead of Eq. (4). Consider only the 
simpler case when /2 = 0. 

Am, H = kp'^x'^WimD''). 
Prob. 13. — Solve Prob. 12 when R is not equal to zero. Hint. — 
In integrating introduce a new variable z = R + xp. 

Ans. H = [K{R + px)'"+i - KR'^^]/(mpD''). 

Problems Involving Two Reservoirs. 

In some cases of water supply more than one storage 
reservoir is provided, and these reservoirs may have to 
be maintained at different levels; for instance, one may 
be the main reservoir and the others auxiliary to it. 
Referring to Fig. 33, let A be the main reservoir which 
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Fig. 33. — A pipe fed from two reservoirs. 



during the hours of light demand suppUes not only the 
entire city but also replenishes the auxiUary reservoir 
jB. During the hours of heavy demand both reservoirs 
discharge into the city mains. This arrangement per- 
mits of a reduction in the size of the pumps, the utiUza- 
tion of the full pumping capacity nearly all the time, 
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and the maintenance of a more uniform pressm-e in the 
mains throughout the day. With this arrangement 
several interesting practical problems arise in regard to 
the interchange of water between the two reservoirs, 
the direction of flow in the connecting pipe, the time 
necessary for a change in the level of the water in the 
reservoirs, etc. 

Prob. 14. — The cross-connecting pipe KL (Fig. 33) is connected 
at C to another pipe in which the demand R is variable. Find the 
relation between this demand and the head h At C; also the parts 
Qi and Q2 of the demand that are supplied from reservoirs A and 
B respectively. The levels at A and B are given and are assumed 
to be constant. Solution. — According to Eq. (4) we have : Ai — A = 
cQiHi/D^, and A2 - ^ = cQ2%/D^\ moreover, Qi + Q2 = -B. These 
three equations contain three unknown quantities, A, Qi, and Qj. 
Eliminating Qi and Q2 we find the following equation for h: 

^{hi-h)/h + y/{h2 -h)/U = v^ . R/D^, 

which is the required relation between R and h. Knowing A, the 
quantities Qi and Q2 are easily determined from the first two equa- 
tions. 

Prob. 15. — In the preceding problem investigate the case when 
h > hiy according to the upper line of charge in Fig. 33. Hird. — 
Reservoir A supplies not only the pipe at C, but also feeds the 
auxiliary tank B, The sign of Q2 in the equations must be reversed. 
In the limiting case when A = fe, the line of charge between B and 
€ is horizontal; this means that no water is supplied from B, and 

Prob. 16. — Apply the formulae of the two preceding problems 
to the case where the difference of the levels between A and B is 
6 m., the pipe is 60 cm. in diameter, 2500 m. long and point C is 
900 m. distant from B. Plot to values of i2 as abscissae, the corre- 
sponding values of loss of head in both parts of the pipe, and the 
rates of discharge Qi and Q2, with the proper sign. On the basis 
of the curves discuss the influence of the auxihary reservoir in as- 
sisting to maintain a better pressure, and in equalizing fluctuations 
in the demand upon the main reservoir. Hird, — The solution de- 
I)ends upon the differences of levels, and not upon their absolute 
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values. The absolute levels limit only the maximum value of R, 
since h cannot be less than — 10.4 m., corresponding to a discharge 
into vacuum. For a discharge into the atmosphere, without any 
free head left, ^ = (see the next problem). 

Prob. 17. — Calculate, as in Prob. 16, the values of Qi and Qa 
when hi = 36 m. and A2 = 30 m.; the discharge R Sii C takes place 
directly into the atmosphere, so that h is nearly equal to zero. 

Ans. Qi = 861 l./sec; Q2 = 1047 l./sec. 

Prob. 18. — The demand, instead of being concentrated at a 
point C (Fig. 34), is distributed uniformly throughout AB, at a rate 
of p liters per meter. Find the discharge Qi and Q2 from both 
reservoirs, the shape of the curve of charge acb and the point C of 




Fig. 34. — A pipe fed from both ends, with a uniformly distributed 
demand. 

division of supply. SoltUion. — Let C be the unknown point of 
division of supply; to the left of this point the demand is supplied 
by reservoir A, to the right by reservoir B, The line of charge has 
a minimum at C, According to Prob. 7, we have for the left- 
hand side hi — h =^ J c (pxY • x/D^y and for the right-hand side 
^2 — h = ic[p (l — x)f(l — x)/D^. Subtracting one equation from 
the other in order to eliminate A, a cubic equation for x is obtained. 
When /t2 = h, the second equation gives j = Z, and the first h — 
^2=1 cpH^/D^, The rate p determined from this equation is 
such that reservoir A supplies the whole demand. If the demand 
falls below this rate, the tank A will begin to supply water to 
tank B. 

' Prob. 19. — Apply the method indicated in the preceding prob- 
lem to the case of a 100 cm. pipe 2000 m. long, the difference of the 
levels between the two reservoirs being 4 m., and the level in the 
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main reservoir 60 m. above the pipe. Find the rates of de- 
mand p at which x = i Z, i Z, and i I; also the minimum rate below 
which reservoir A not only supplies the whole demand but also 
replenishes tank B. 

Prob. 20. — The demand is imiformly distributed (Fig. 34), as in 
Prob. 18, and in addition there is a concentrated demand R &i & 
given point F of the pipe, determined by the distance m. Find an 
expression for the position of the point C of division of the load, and 
the equation of the line of charge. 

Am, Um>x; h-hi^ {c/D^) {iT^-phnx^+P^'h^-M 

+R^-Rpl^+2Rplm-Rpm^-j}^P/3+j)^Pm-'p^W) ; 

Um<x; hi-h2= {c/D^)(j)h^+Rpx^-p^lx^+R^x+p^i:'x 

- p^yS - p2p/3). 

Prob. 21. — The original difference of levels in the reservoirs 
(Fig. 33) is h\ Find the difference of the levels h" after an interval 
of time T, when the water is allowed to equalize freely between the 
two reservoirs and there is no demand in the branch pipe at C. 
Solution, — Eq. (4) gives 

h^hi-h^ m/J^\ (10) 

where hxy h, and Q are functions of time, so that the equation 
applies only to an infinitesimal interval of time dt. As much water 
flows out of A as flows into B, so that 

-Fidhi = F2dh2 = Qdt (11) 

In these equations Fi and F2 stand for the constant areas of hori- 
zontal cross sections of the reservoirs. The minus sign before dhi 
is necessary because hi decreases with time so that dhi is a negative 
quantity; the minus sign makes it positive. Otherwise it could 
not be equal to F2 dh^ and to Q eft, both of which are essentially posi- 
tive quantities. From Eqs. (10) and (11) the variables Ai, ^2, 
and Q may be eliminated and the resultant equation containing h 
and t can be integrated. 

Differentiating the first two members of Eq. (10) and comparing 
them with the first two members of Eq. (11), we get: dhi — FidhJ 
(Fi + F2), while from Eq. (10) Q = D* . >/h/y/d. Substituting 
these values in Eq. (11) we obtain 



Fi + F, Va 

/Google 



Digitized by ^ 



70 ENGINEERING MATHEMATICS. [Chap. V. 

We introduce, for the sake of abbreviation, F = FiF^/^Fi + Fj), 
and a = D^/Vd, Then 

- f^''dh/Vh= {a/F) Cdt, 

or finally, 

VA'- VF = (a/2iJ').r (12) 

The final difference of levels h" is determined from this equation. 
If h!* = 0, Eq. (12) gives the amount of time necessary for equaliz- 
ing the levels in the two reservoirs. 

Prob. 22. — In the preceding problem determine the changes in 
the levels, Mi and AA2, in the two reservoirs. 

Am, AAi = F2 {K - h!')l{F^ + Fi), 

Prob. 23. — Solve Prob. 21, using Eq. (5) instead of Eq. (4). 

Am. Qi'Y- {h'y=-qbT, whereg= l-m-^, and 6 = DV"»/(A;0i/« 

Prob. 24. — Solve Prob. 21 when a quantity of water Qi per sec- 
ond is pumped into reservoir A, and a quantity Q2 is withdrawn 
from reservoir B. Hint. — Using Eq. (11) we obtain: —iJ'i d/ti 
= (Q— Qi) dtf and F2dh2 = (0 + Q2) dt. After eliminating hi and 
A2 the result is 

j,^_ r^" FiF2dh 

Jh^ a (Fi + F2) Vh + {F1Q2 - F2Q1) ' 
where a = D^/Vd. The integral is reduced to a simple standard 
form by introducing a new variable z = Vh. 

Pipe Line of Vaiying Diameter. 
Prob. 25. — A street main has been gradually extended, so that 
at a certain time it consists of sections of different lengths h, fe, etc., 
and of diameters Di, D2, etc. It is desired to replace this by a new 
pipe of a uniform diameter that would give the same loss of head at 
the end, with the same rate of discharge. Find the diameter Z)eq 
of this equivalent new main. Solution. — According to Eq. (4), 

from this relation Deq may be calculated. If Eq. (5) be used, the 
exponent 5 must be replaced by n. 

Prob. 26. — A conical pipe has a diameter A at one end and D% 
at the other end. The length of the pipe is I. Find the expression 
for the loss of head. Solution. — According to Eq. (4) 
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where the diameter at the distance x from the end Di is 

D = Z)i + (A - DO .x/Z, (13a) 

80 that dD = (Z>2 — Di) • dx/l. Substituting, we get 

D2-D,Jd, D^ 4 (A -A) I a* a*)* • ^ ^ 

This formula is inconvenient for computations when D2 differs but 
Httle from Di, because the differences (D2 — Di) and (D2* — A*) 
are very small, and a minor inaccuracy in each may affect the result 
to a considerable degree. In this case put A = A (1 + a), where 
a is a small quantity as compared with 1. Substituting we obtain 

4A^-a( (1 + a)^ 

Expand (1 + a)"* according to the binomial theorem, subtract the 
series from 1, and divide the result by a. Then 

A = ^^jl_2.5a-5a2- . . .}. . . . (15) 
JJi° 

The correction a being a small quantity, the rest of the terms of the 
series may be neglected. The expression in the brackets shows the 
effect of tapering the pipe, as compared with a cylindrical pipe. 

This method of transformation will be found useful in other 
problems in physics and in engineering, when one has to subtract 
two quantities nearly equal to each other. 

Prob. 27. — Investigate the relative advantages and conven- 
ience of formulae (15) and (14) when A is greater than A by 1, 3, 
10, and 100 per cent respectively. 

Prob. 28. — A conical pipe 100 m. long has a diameter of 5 cm. 
at one end and 12 cm. at the other end. What is the diameter of 
an equivalent cylindrical pipe, that is, of a pipe of the same length 
which gives the same total loss of head at the same rate of dis- 
charge? Ans, 7.1 cm. 

Prob. 29.' — Deduce the equation of the line of charge for a con- 

1 1 1 



icalpipe. . , cQH 



A* /n , A - A 



(z>,+5^.«) 



Prob. 30. — Plot the line of charge for the conical pipe given in 
Prob. 28, and for the equivalent cylindrical pipe, the total loss of 
charge being in both cases 0.2 m. 
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Prob. 31, — For a conical pipe deduce an expression correspond- 
ing to formula (9). Hird. H = cf[R + p{l- x)Y dx/D\ 

•'o 

where D is given in Eq. (13a). In order to perform the integration, 
take D for the independent variable, and express z in terms of Z). 

Prob. 32. — Show how to modify formulae (14) and (15) when the 
increase in the diameter of the pipe is proportional to the square of 
the length, so that Eq. (13a) becomes 

Z) = Z)i+(Z>2-A)(x/Z)2 (16) 

Also extend the treatment to the case when the increase in the diam- 
eter is proportional to some nth power of the length, or in general 
when D is a given function of x. 

Problems Involving Minimum Expense. 

In deciding upon the sizes of pipes and pumps, and 
upon the general layout of a water supply system, the 
commercial question of minimum expense enters as 
one of the most important factors. To design a water 
supply system, including the pumping station, the reser- 
voirs, and the mains, so as to have the most economical 
installation, both as to first cost and cost of operation, 
is the problem that confronts the engineer and at the 
same time makes his task much more definite than it 
would be otherwise. The economy problem is different 
in different cases; for instance, it may be required to 
design a plant in such a way as to obtain minimum first 
cost of installation irrespective of operating expenses; 
or minimmn operating expenses may be desired irre- 
spective of first cost; sometimes maximiun net revenue 
is important, etc. 

It is hardly possible or even desirable to express in 
the form of equations the practical problem in all its 
complexity. Moreover, the values that correspond to 
an absolute minimum, according to the rules of the dif- 
ferential calculus, may not be the best, considering the 
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fact that the engmeer is Umited to comparatively few 
standard sizes of pipes and pumps. However, the vari- 
ous parts of the problem may be worked out mathemati- 
cally and thus furnish at least a first approximation as to 
the size of the pumping engines, diameter of the mains, 
etc. The method is explained below with the aid of a 
number of problems graded from the simplest to those 
comparatively compUcated. Many of these problems 
are appUcable with certain modifications to electrical 
distribution; conductors being substituted for pipes, 
and generating sets for pmnping engines. See Chapters 
II and III, in Part V of this work. 

The fundamental assumption made is that the imit 
prices of machinery and pipes are known functions of 
their sizes. Thus, the cost of an installed pipe includes 
the cost of the pipe itself, the work of making the 
trench and filling it, of connecting the sections, trans- 
porting the pipe, etc. Part of this expense, such as the 
cost of the pipe itself and its transportation, depends 
upon the size of the pipe; other items of expense, for 
instance digging the trench, are practically independent 
of the size of pipe. Thus, in the simplest case the cost 
per meter of installed pipe may be expressed by a + bD, 
where a represents the expense independent of the size 
of the pipe, and bD is the part proportional to the 
selected diameter D of the pipe. The cost of installa- 
tion of a system of pipes is therefore Si (a + 6Z)), and 
this expression must be a minimum; or 

dEl + hZlD = min. 

In most cases the lengths or at least the smn of the 
lengths of the pipes is a fixed quantity, so that the 
above condition is reduced simply to 

JllD = min (17) 
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Prob. 33. — A pipe line of length I must supply a given quantity 
of water Q per second, at a given loss of pressure h. It is proposed 
to make it out of two sections in series, so as to reduce the expense 
to a minimum. Find the diameter and the length of each section. 
SoliUion. — In order to obtain a more general answer Eq. (5) will 
be used. Let Zi, k, Di, and D2 be the imknown lengths and the 
diameters of the sections. The conditions of the problem may be 
expressed thus: 

^=^«1^+^! (^^> 

h + k^l (19) 

h (a + bDi) + k{a + hD^) = min. 

By omitting the constant part, the last equation may be simplified 
to an equation corresponding to Eq. (17), thus: 

ZiZ)i + y)2 = mm (20) 

The answer to the problem may be foreseen directly from the fact 
that the two variables Di and Z>2 enter into the equations symmet- 
rically, and so do the other two variables h and Z2. The answer 
therefore must be Di = A and ^1 = ^2; in other words, pipe of one 
diameter is cheaper than that of two, whatever their diameters 
may be. 

If, however, an exact solution is required, it may be obtained 
according to the general rule for finding relative maxima and minima. 

Multiply Eqs. (18) and (19) by indeterminate coeflScients Xi and 
X2, and add to Eq. (20). The result is 

Partial derivatives with respect to A, A, ^i, and fe, equated to zero, 
give: 

-^. + '^ = «' ;^ + x, + A = 0; 

-^ + ^ = «' ^„ + x, + A = 0. 

From the first two equations, \x/D{^ = Di/n, and Xi/A" = A/n. 
Substituting these values in the last two equations, we easily ob- 
tain Z)i = D2,or the pipe must be of the same diameter throughout. 
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A similar treatment applies to electrical conductors in series; see 
Part V, Chap. II, Prob. 1. 

Prob. 34. — A 60 cm. pipe supplies water to a reservoir, with a 
given loss of head. What would be the extra cost of the pipe line if 
one half of the pipe were replaced by a larger pipe Di cm. in diameter 
and the other half by a smaller pipe so as to have the same total loss 
of head as before? Plot a curve showing extra expense for Di equal 
to 60, 70, 80, 90, and 100 cm., assimaing the expense for 60 cm. pipe 
to be 100 per cent; also plot the corresponding values of A. Con- 
sider the total expense to be represented by formula (20) ; use for- 
mula, (4) instead of formula (5) . Hint. (60)-« = 0.5 (Dr* + Di"') ; 
per cent extra expense = 100 (0.5 A + 0.5 A - 60) /60. 

Prob. 35. — Two pipes in parallel, of diameters Di and D2 re- 
spectively, supply a reservoir. Show that for the same service one 
larger pipe would be cheaper. 

Solution. — Let the diameter of the equivalent pipe be Deq. 
Then we have, according to Eq. (4): Qi = VDi^ • Vtjc; Q2 = 
vW- v^; Q = VD^^ . V?7^, and Q = Qi + Q2. 

From these equations we obtain by adding the first two and 
equating their sum to the third: 

V^ = Vd? + VD^. (21) 

This equation determines the diameter of the equivalent pipe. 
In order to prove that, one pipe is cheaper than the two, it is neces- 
sary to demonstrate that 

2>eq<A + A (22) 

To do this, raise both sides of Eq. (21) to the power of 2/5, using 
the binomial expression. The result is that Deq = Di + Z>2 + other 
positive terms, whence inequality (22) follows directly. 

It there are more than two pipes in parallel, the cost of installar 
tion may be reduced by replacing any two by a larger pipe; then 
again this pipe may be combined with another, and so forth, so 
that fi nally one pipe, whose diameter is determined from the rela- 
tion VDeq* = SVZ)<>, is more economical than any given system 
of pipes in parallel. A similar result may be proved for electrical 
conductors in parallel. 

Prob. 36. — A 60 cm. pipe supplies a reservoir. Show how much 
more expensive the installation would be if there were two supply 
pipes in parallel, of which one might be of diameter A = 50, 40, 
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30, 20, 10 cm., and the other of such a diameter D2 as to give the 
required discharge at the specified loss of head. Assume the cost 
of the original pipe to be 100 per cent, and represent the results by 
curves. 

Arts, For A = 30 cm., A = 55.5 cm., cost = 142.5 per cent. 
Prob. 37. — A water main carries Qi liters per sec. for a distance 
li (Fig. 35), and Q2 for a distance U- Determine the diameters of 



! r> i n j 



Dj \ D, 



~. I =^' 1 • 

Fig. 35. — A composite water main \vith a branch pipe. 

the two parts of the main so that its cost will be a minimima, for a 
given total loss of head. Solution. — 

A^ ^ W ~ ' ^ ^ 

ZiDi + Z2A = mm (24) 

The complete differential of the left-hand member of the first 
equation is equal to zero, because A is a constant; the complete dififer- 
cntial of the left-hand member of the second equation is also equal 
to zero, this being the condition for a minimum. Thus 

liQi^ dDi/Di^ + m2^ dD2W = 0, 
and 

lidDi + UdDz^O. 

Substituting h dDi from the second equation into the first we get 
after reduction: 

Di/D2=\^Q^2; (25) 

this equation, together with Eq. (23), allows one to determine the 
unknown diameters A and A- 

Prob. 38. — Show that in the preceding problem the hydraulic 
gradient is greater in the first pipe than in the second, and that for 
a minimum cost the diameters of the pipes must be as their hy- 
draulic gradients. 
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Prob. 39. — Solve Prob. 37 using Eq. (5) instead of Eq. (4). 

Prob. 40. — Extend Prob. 37 to the case when there are more 
than two sections in series. Solviion. — 



^XS = ^' (26) 



SZZ) = min (27) 

According to the general rule for relative maxima and minima, we 
form the expression 

" = (21-!)+^"'. 

where X is an indeterminate coefficient. The partial derivatives 
of B with respect to all D's must be separately equal to zero. All 
these derivative equations are the same form: 



^+XZ» = 0, 



from which 



^' = i/5 = const., 



so that 



0- 

This expression is analogous to Eq. (25). Coefficient X is determined 
by substituting the values of D's from Eqs. (28) in Eq. (26). 



|-« 50 111. »^<-20Tn.->|< — 40 m. — >*< 60 m-, >' 

1 ' " ; ^ _i 

loysec. ^Veec, 15^/sec. 8?/sec. 

Fig. 36. — A street main of tapering diameter. 

Prob. 41. — Determine the most economical sizes of pipe for 
the main shown in Fig. 36, when the total permissible loss of head 
is 10 m. The amounts of discharge shown in the figure represent 
deUvery to side streets. Thus Qi = 10 + 4 + 15 + 8 l./sec; 
62 ^ 4 + 15 + 8 l./sec; etc. 

Ans, 2)i = 12.54cm.;2)2=11.91cm.;Z)j=11.6cm.;Z)4=9.72cm. 
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Prob. 42. — Calculate the three diameters of a branched mam 
(Fig. 37) to correspond to a minimum expense of installation. The 
rates of discharge and the total permissible loss of head h from A to 
B and from A to C are given. Solution. — Let the unknown loss 
of head from A to iiT be x. Then we have four conditions: x = 



I 



Qi + Qs 




Fig. 37. — A branched main. 

« (Oi + Q2yi/D^; (A - a:) = cQiHi/Di^; (h - x) = cQz^/W; and 
ID + liDi + U)2 = min. Substituting the values of D, Di, and 
D2, from the first three equations in the last, we obtain 

^^ 7 1 7 H 7 = mm. 

a:* Qi - x)* Qi - x)^ 

Equating the first derivative to zero, and rearranging the terms we 
get 

x* Qi - x)^ Qi - x)* 

from which x = h/(l + n), where 

^^lOM + Qll^. 
l(Qi + Q2)^ 

Prob. 43. — Extend the result obtained in the preceding problem 
to the case when there are several branches at K, instead of two; 
also when formula (5) is used instead of formula (4). 

Arts. X = h/(l + n) where n = — - — 7- • 

Prob. 44. — Referring to Probs. 42 and 43, and to Fig. 37, let 
the total loss of head be different for each point, such as B, C, etc. 
Show that the problem can be solved only by trial and successive 
approximations, and indicate a convenient order of computations. 
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Prob. 45. — Solve Prob. 43 using the method of indeterminate 
coefficients, as in Prob. 33, and show that it leads to the relationship 

Prob. 46. — Establish equations for minimum expense for the 
branched pipe line shown in Fig. 38. Hint, — As the independent 




Fig. 38. — A branched pipe line. 



variables, select the loss of head x from A to K, and the loss of head 
y from K to L. 



+ 



and 



y {he -X- 2/)* QiD-x- 2/)*' 

2/* X* (hs — x)^ 



Prob. 47. — Referring to the preceding problem, find the pipe 
diameters for the following conditions: Qb = 251. /sec; Qc = 10 
l./sec; Qd = 60 l./sec; I = 50 m.; Zi = 70 m.; Z2 = 160 m.; 
k = 120 m.; ^4 = 70 m.; hs = 10 m.; he = ho = 12.5 m. 

Am. D = 22.1 cm.; A = 18.1 cm.; D2 = 7.8 cm.; Ds = 8.5 cm.; 
Di = 15.7 cm. 

Prob. 48. — A reservoir (Fig. 39) is to be fed by gravity from the 
river AB through a pipe AM, Find point A at which the pipe must 
originate in order to make its cost a minimum. By making the 
pipe longer, the hydraulic gradient i is increased so that the diameter 
of the pipe may be reduced for the same Q, Since the cost of the 
pipe is proportional to DZ, there must be values of D and I which 
will make this expression a minimum. Solviion, — Let I be the 
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slope of the river, or the trigonometric tangent of the angle which 
its surface makes with the horizontal. Then we have from the 
triangle AMN: 

A = Z(7-i) (29) 

Besides 

i^dQ'/D^ (30) 

and 

ZD = min (31) 

Eliminating I and D from Eq. (31) by means of Eqs. (29) and (30) 
we find that i~* (/ — i)~^ must be a TniniTnnm , or i* (7 — i) must be 




Fig. 39. — A reservoir fed by gravity from a river. 

a maximum. According to the familiar rule of the differential 
calculus we find that i = J 7. This determines the length and the 
diameter of the pipe. 

Prob. 49. — In Prob. 48 the cost of the trench and condemnation 
of the ground is to be taken into account. Solution, — Let p be 
the cost of the ground and excavation per linear foot of the pipe 
line; this cost is practically independent of the diameter of the 
pipe. Let q be the cost of the pipe per linear foot and per unit 
diameter. Then expression (31) becomes 



{qD + p) Z = min. 



(32) 



Eliminating I and D as before, an equation is obtained for i which 
is then solved for a minimum. In this case it is simpler, however, 
to eliminate i and Z, and to establish an equation with respect to 
Z). Before doing this, it is convenient also to eliminate Q, and thus 
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to make the equation more symmetrical by introducing a known 
fictitious diameter A such that 

/ = cQ2/A« (33) 

Eq. (30) becomes 

t = I(A/Dy (34) 

Substituting the values of I and i from Eqs. (29) and (34) into Eq. 
(32) we obtain, after omitting the constant factors: 

* = mm. 



-(?)■ 



In this equation it is convenient to consider the ratio (D/A) aa the 
variable; differentiating with respect to this variable and equating 
the derivative to zero gives the following equation of the sixth 
degree: 

(!)'-H!)-^-« '-' 

This equation can be solved graphically, or by the method of succes- 
sive approximations. 

Prob. 50. — Show that in Eq. (35) the value of (D/A) is always 
larger than 1.43. Hint, -^ Put p = 0, 

Prob. 51. — Using Eq. (35), plot a graph of values of (D/A) 
corresponding'to values of p/qA from to 2. Analyze the practical 
significance of such a curve. 

Prob. 52. — A pmnping station must supply Q liters of water 
per second to a reservoir at a height H above a river and distant I 
meters from it. Determine the capacity of pumping engines and 
the diameter of the pipe such as to reduce the first cost of the whole 
installation to a minimum. Solviion, — The output of the pump- 
ing engines is spent partly in lifting the quantity Q of water to 
a given height H and partly in overcoming friction in the pipe. 
The first item is fixed, the second depends upon the diameter of the 
pipe. If a smaller and thus less expensive pipe be selected, more 
work will be required for overcoming friction, and vice versa. The 
loss of head, according to formula (4), is A = cQ^l/D^, and the power 
wasted, corresponding to this loss, is cQ^l/(D^ri) kg. met. per second, 
where 17 is the mechanical efficiency of the engines. 
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Let the cost of one additional indicated unit of power of the 
engines, boilers, pumps, foundations, etc., be p dollars, and let the 
cost of pipe per linear meter and per centimeter of diameter be q 
dollars. Then the condition of minimum cost becomes 

pcQ^l/iD^v) + qlD = min. 

This expression becomes a minimum when D = VQ . -^5 pc/iqri). 
The student should satisfy himself as to the physical reason for 
which this value of D is independent of I, 

Prob. 53. — In the preceding problem determine the size of the 
pipe in such a way that the total operating expense, including the 
interest and depreciation on the plant, will be a minimum. Solvtion, 
— If a be the rate of interest and depreciation for the pumpingnsta- 
tion equipment, h that for the pipe, and / the cost of operating the 
plant per extra unit of power per year, the expression for the operat- 
ing expense becomes 

cQ^l (flP +f)/D^v + hqlD =;min., 

which is solved for D as before. 

Prob. 54. — How would the foregoing expression be modified 
if the plant were operated only one nth part of the year, instead 
of the whole year? Would a larger or a smaller diameter of pipe be 
selected, as compared with the pipe selected when the plant is in 
continuous operation? 
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CHAPTER VI. 

BEST FORM OF CHANNEL SECTION. 

When an open channel, such as an aqueduct or a 
sewer, is arranged for a flow of water, the bottom of 
the channel must be given a longitudinal slope in order 
to overcome the friction of water against the walls. 
It is desirable to give the cross section of the channel 
such a form as to reduce this friction to a minimum. 
The quantity of water Q discharged per minute is usually 
given, as is also the desirable average velocity v of flow. 
Hence, according to the equation 

Q = vA, (1) 

the area A of the cross section of the channel is a definite 
given quantity. The problem is to find such a shape 
of the cross section that the friction will be a minimum. 
The friction of liquids against a solid wall is propor- 
tional to the wetted surface, or, per unit length of the 
channel, is proportional to the wetted perimeter of 
the cross section. Hence, to reduce the friction, the 
wetted perimeter must be kept as low as possible. 
The problem is thus reduced to the following geometrical 
proposition : 

To find the form of a cross section which for a given area 
possesses a minimum perimeter. 

Of all possible geometric forms the circle has the 
absolute minimum of perimeter for a given area. For 
this reason a round pipe is the most rational form 
where it can be used. For an open canal, a semicircular 
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cross section is just as good, because while the area 
is halved the wetted perimeter is also cut in two. 

It is not practicable, however, on account of expense, 
to provide canals which are dug out in the ground, 
or even those lined with masonry or concrete, with a 
semicircular cross section. Next to the circle, regular 
polygons (Fig. 40) possess a comparatively small perim- 
eter for a given area, the perimeter being the smaller 
the larger the number of sides, in other words, the more 
closely the polygon approaches the circle. Thus, when 



Fig. 40. — Some cross sections of channels. 

a pipe cannot be made round, for instance when it is 
made out of boards, like big penstocks for hydraulic 
plants, it should be made hexagonal or square. In the 
case of open channels, the corresponding cross sections 
are found by taking the lower part of the figures (Fig. 
40), because the area and the wetted perimeter are 
reduced in the same proportion. 

There are many practical cases in which regular cross 
sections cannot be used, for instance, when the slope 
of the sides of a channel is prescribed in order to keep 
them from sUding, when the depth of a canal is pre- 
scribed for navigation, etc. In such cases the shape 
of the cross section may be found according to the 
general rule given in the differential calculus for finding 
relative maxima and minima. 
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Prob. I. — Prove that the best form of a rectangular cross 
section for an open channel is one with the width equal to twice 
the depth. 

First SolvJtion. — Let the depth of water be y\ then the width 
of the channel is A/t/, where A is the given area of cross section. 
Hence the wetted perimeter P = 2y + A/t/ = min. This ex- 
pression becomes a minimum, when dP/dy = 2 — A/y^ = 0, or 
y = y/iA. The width A/y is equal to V2 A or 2 1/. 

Second Solviion. — Let the width be x and the depth y] then we 
have two equations : 

A = xy = const.; 

P = 2t/ + x = min.; 

dA = 0, because A is const. ; dP = 0, because P must be a mini- 
mum; so that differentiating the two equations we get 

xdy + ydx = 0] 
2dy + dx = 0. 

Eliminating dx and dy from these equations gives x = 2 y, as 
before. 

Theorem 1 in the Appendix may be applied in this case, because if 
xy = const., 2y •x is also a constant. 

When for practical reasons it is not possible to give a 
channel the theoretically best possible cross section, it is 
desirable to ascertain at least in how far the shape actu- 
ally selected departs from the one most advantageous. 
This means, how much greater longitudinal slope 
must be given the channel in order to obtain the same 
velocity of water as with the theoretically best form of 
cross section, or how much larger the cross section must 
be for the same discharge with the former slope, or how 
much less water will flow with the same area of cross 
section. The relationship between the quantities in 
the problem is obtained as follows: 

Imagine a prism of water of length I and of cross 
section A (Fig. 41), which slides along a canal. The 
motion being uniform, the component of the weight 
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of the water along the bed must be equal to the retard- 
ing force of friction. We have thus 

Alw sin a = fPL 

In this equation Alw is the weight of the water, a 
is the angle of inclination of the bed to the horizontal, 




Fig. 41. — Motion of water in a canal. 

/ the force of friction per unit area, and PI is the wetted 
surface. 

Experiment shows that the unit force of friction / is 
nearly proportional to the square of velocity, or 

/ = const. X v^. 

Angle a is usually so small that instead of sin a we may 
substitute tan a = i, where i is the slope of the bed. 
We thus obtain 

Awi = const. X v^P, 

or, collecting the constants into one constant c, we 
finally get 

v^cViAjP, (2) 

which is the fundamental formula for the flow of water 
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in channels.* Combining Eq. (2) with Eq. (1), another 
fundamental formula is obtained. 

Q = cAViA/P (3) 

Prob. 2. — A rectangular canal, in order to be of the most ad- 
vantageous form, should have a depth h and a width 2h, For 
certain practical reasons the depth has to be changed to ph, where 
p is a certain coefficient which may be greater or smaller than unity. 
It is desired to preserve the same cross section 2 h^ and the same 
rate of flow of water, so that the longitudinal slope i must be in- 
creased q times. Find the relationship between p and g, plot the 
curve between them and define the curve mathematically. Solu- 
tion. — If the new depth is ph, the width must be 2 h/p, and the 
new wetted perimeter P^ = 2ph + 2 h/p. But, according to for- 
mula (3), both the new and the old values of the slope must satisfy 
the relation V/i = P^Pj or 

^ 2ph + 2 h/p 
^ 2h + 2h 

Simplifying, this gives the required relationship 



-W-i)- 



If, for instance, the depth had been increased by 20 per cent, we 
would have 



*=l(^-2+^)=i-«i«' 



or the new slope must be 1.6 per cent greater than the old. It is 
left to the student to plot the curve and to prove that it is an hyper- 
bola; also to find the direction and the magnitude of its axes. 

Prob. 3. — A channel should have a depth h and a width 2 h. 
For certain reasons the depth must be changed to phj and in order 
to preserve the same rate of discharge at the same slope the area 
of cross section must be increased r times. Find the relationship 
between p and r. Ans, r* = i (p + r/p). 

Prob. 4. — If in the preceding problem the area be kept constant, 
the discharge decreases, so that if the original rate of discharge is 
g, the new rate iaq' = itq. Find the relationship between p and u. 

Ans, ^2 = 2p/(p2+ 1). 

* In reality c varies somewhat with many factors, but in these 
problems it is assvuned to be constant. 
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Prob. $. — In a rectangular canal of width b, the depth of water 
is different at different times. Deduce an equation and plot a curve 
between the depth of water h and the average velocity of flow v. 

Assume the velocity at A = i 6 to be unity. Ans, t^ = ^ , ,,, ' 

2 n/o -f- 1 

Prob. 6. — In a rectangular canal of width 6, the depth of water 

is different at different times. When h= ih,Q = Qo, Deduce the 

equation and indicate the general character of the curve between the 

depth of water h and the corresponding rate of discharge Q. 

Ans -0!__16W^. 
^'^' Qo'~2h/h + l 

When a canal is dug in the ground, it usually has a 
trapezoidal cross section (Fig. 42) in order to prevent 





Fig. 42. — A trapezoidal canal. 

its banks from sUding. The area of the cross section is 
A = y {x +y tan e), .... (4) 
and the wetted perimeter is 

P=x+-^ (5) 

cos e 

The cross section is determined by three parameters, 
X, y, and 0, and as long as none of them are given, the 
problem of designing the best cross section remains 
indefinite. The two conditions, A = const, and P = 
min., are not sufficient to make the problem definite. 
In practice, however, either 6 is given, according to 
the character of the soil, or of the Uning, or else the 
depth y is prescribed. Sometimes the minimum width 
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on the bottom is prescribed, for instance, of such a size 
that two barges may pass by one another. We shall 
consider these cases separately in some of the following 
problems. 

Prob. 7. — Find the best dimensions for a trapezoidal canal, the 
slope of the banks being given. Solviion. — In Eqs. (4) and (5) 
angle e must be considered constant, and we have: 

dA = {x + y tan e)dy + y{dx + tan e^dy) = 0; 
dP = dx + 2dy/ cos ^ = 0. 

Eliminating dx and dy from these equations gives one equation 
between x and y. Relation (4) is the other necessary equation. 

Am. y^ = A cos e/{2 — sin ^). 

Prob. 8. — On the same base line and to the same scale draw 
trapezoidal cross sections of canals for minimum friction and for 
the slopes of the banks equal to 2 : 1, 1.5 : 1, 1 : 1 and 0.5 : 1, respec- 
tively. The area of the canal is to be the same in all cases; use 
the answer to the preceding problem. A slope of 2 : 1 means that 
for every two feet of horizontal projection the surface rises one 
foot. 

Prob. 9. — Solve Prob. 7 when the width at the bottom is given 
instead of the slope. Solution. — DiflFerentiating Eqs. (4) and (5), 
and considering x to be constant, we obtain: 

dA = {x + y tan e)dy + y (tan edy + y dd/cos^ 6) = 0; 
dP = 2 dy /cos ^ + 2 y sin ^ de/cos^ 6 = 0. 

Eliminating dy and de from these equations gives an equation be- 
tween y and By which being solved with Eq. (4) gives the values of 
the imknown quantities. It is left to the student to obtain the 
final expressions f or y and 6. Ans. y^ + Axy = A^. 

Prob. 10. — Solve Prob. 7 when the depth y of the channel is 
fixed, instead of the slope. Ans. x = A/y — y/VS, 

Prob. II. — Express Eqs. (4) and (5) directly through the slope 
n = tan ^, and solve Probs. 7 to 10 in this form. 

Prob. 12. — Solve Prob. 7, when the slope of the two banks is 
different. Ans. 2 A = z/^ } (2 — sin ^i) /cos ^1 + (2 — sin ^2) /cos $2 \ . 

Prob. 13. — For a trapezoidal channel with a given width at 
the bottom and with a given slope of the sides (the slope may be 
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different for the two banks) find the depth of water at which the 
average velocity of flow is a maximum. 

2xy 2x2 



Arts, y^ + " 



= 0. 



sec $1 + sec 62 (tan $1 + tan $2) (sec ^1 + sec ^2) 

Prob. 14. — Solve Prob. 13 when the bottom line of the cross 

section is inclined at a given angle a (Fig. 43) to the horizontal. 

A •10 1 — sin a sec ^2 

Ana. y^ + 2xy z ; ^— r 

cos a (sec ^1 + sec ^2) 

, , (1 — tan a tan ^2) (2 +sin a sec gi — sin « sec $2) _ ^ 

cos a (tan $1 + tan $2) (sec Si + sec ^2) 




Fig. 43. — A canal in which the bottom line of the cross section 
is inclined at an angle a to the horizontal. 

Prob. 1$. — A certain quantity of water is conducted away by a 
trapezoidal canal of minimum frictional resistance. What would 
be the per cent saving in excavation if the channel could be given a 





Fig. 44. — A partly filled circular Fig. 45. — A sewer with a square 
sewer. cross section. 

semicircular cross section? The longitudinal slope is supposed to be 



the same in both cases. 



Ans, 100 



ix + 2yBece) 



100. 



• 2 y sec d) 

Prob. 16. — For a partly filled circular sewer (Fig. 44) find the 
level of water at which the velocity is a maximum. SoltUion. — 
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Let the level be determined by the angle 0, Then the wetted area 
A = area of the circle — area of the upper segment = irr* — i r* 
(^ — sin 0). The wetted perimeter P = r (2 x — ^). The ratio of 
A to P must be a maximum, or d (A/P) = 0. This leads to the 
condition A dP - P dA = 0. But rfA = -i r* (1 - cos d) d0, and 
dP = — r d0, so that 

[irr*- ij^ {0 - em0)]rde ^ r (2t - 0)'ir'^ (1 - co8^)d^, 

whence sin ^ = — (2 x — ^) cos d, or tan = — (2 x — ^). The angle 
which satisfies this condition is 102.5°. 

Prob. 17. — Assimie the average velocity of water corresponding 
to the level in Prob. 16 to be 100 per cent, and plot a curve of 
average velocities af other levels of water. 

Ana, ^ = 270 180 90 degrees; 
t;= 62.1 91.1 98.1 91.1 per cent. 
Prob. 18. — Solve Prob. 16 when the level must be such as to give 
a maximum rate of discharge. Hint. — According to Eq. (3) , A'/P 
must be a maximum. Ana. — 52®. 

Prob. 19. — Assuming the maximum discharge deduced in the 
preceding problem to be 100 per cent, plot a curve of rate of dis- 
charge at other levels. 

Ana. ^ = 270 180 90 60 degrees; 
Q= 5.21 47.8 95 99.5 95.3 per cent. 
Prob. 30. — Solve Probs. 17 to 20 in application to a square 
cross section with a vertical diagonal (Fig. 45). Hint. A — a^ — 
jT*; P = 4a-2a: V2(froma: = a/v^ tox = 0). 




Fig. 46. — A circular aqueduct and a corresponding rectangular canal. 

Prob. 31. — A circular aqueduct (Fig. 46) of diameter'rf = 120 
cm. carries water to a level h = 100 cm. For a part of the length it 
is necessary to change the section to an open rectangular section. 
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Calculate the width x of this section such that at the same slope the 
water would flow at the same height h. Hint. — The expression 
A^/P must be the same for both cross sections. Ans. 104 cm. 

Prob. 22. — An egg-shaped sewer (Fig. 47) is determined by the 
dimensions r = 0.25 m., R = 0.585 m., and Z = 1.5 m. The 
middle part is formed by two arcs of circles tangent to the upper 
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Fig. 47. — An egg-shaped sewer. 

circle at the horizontal diameter and to the lower circle. Find 
algebraic expressions for, and plot curves of A and P f or different 
levels of water above and below the center of the upper circle. Find 
whether the ratio A/P, and consequently the average velocity of 
flow, remains sensibly constant within wide limits of water level. 
Ans. H = 1.0; A = 0.38; P = 1.22; A/P = 0.312. 
H = 2.0; A = 0.85; P = 2.46; A/P =^ 0.345. 
Prob. 23. — The cross section of a canal is a cycloid (Fig. 48) 
formed by rolling a circle of radius p on the base line MN. The 
equation of the cycloid is x = p (a — sin a), and t/ = p (1 — cos a), 
where a is the variable angle of rotation of the circle. Find the 
area A of the cross section MNTy the wetted perimeter P, and the 
ratio A/P. Am, A = 3 irp*; P = 8 p; A/P = 3 irp/8. 
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Prob. 24. — Find analjrtic expressions for, and draw curves of 
velocity and discharge at different partial levels of water for the 
cycloidal section in the preceding problem; assume the velocity 
and the discharge at full section MN to be 100 per cent. Hint. — 




Fig. 48. — Cycloidal section of a canal. 



The area of an infinitesimal horizontal strip is 2 (irp — x) dy; the 
limits of integration are y and 2 p, or a and tt. In integra ting to 
find the perimeter rem ember tha t sin i a = V^ (1 — co s a), 
Ans. t; = 46.1 Vm/cos J a; Q = 9.76 Vm'/cos i a, where 
m = a cos a — H*" — a)— X cos a — sina— isin2 a. 




Fig. 49. — Section of a canal formed by a semicircle and two catenaries. 



Prob. 25. — The section of a canal is formed by a semicircle of 
radius 2 a (Fig. 49), and by parts of catenaries MRN and M'R^N\ 
whose equation is 

where e is the base of the natural logarithms. Show that as long as 
the level of water remains above the center of the semicircle, the 
mean velocity of flow is constant. 
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Solution. — The mean velocity depends on the ratio of A to P, 
For the catenary part of the section we have: 

dAc = 2xdy = a {e^l'' + e-^l'') dy, 

dPc = 2 V(dyy -f (dxy = 2 V{dyy + i (e^/" - e-^/y {dyY 
= 2 Vl + i(e2«'/« + e-^''/«-2) .dj^, 

or dPc = (e^/* + «-«'/'•) d2/. 

Consequently, dAc = a rfPc, or, integrating between the levels LU 
and MM'y we find that Ac — aPe- Therefore, for the whole cross 
section, 

iVofe 1. — The problem is more readily solved by using hyper- 
bolic fimctions. The equation of the catenaries becomes ± x = 
a cosh (y/a)j and we have: 

dAc = 2 a cosh (y/a) dy; 



dPc = 2 Vl + {dx/dyy^dx = 2 Vl + smh^ (y/a) .dx= 2 cosh (y/a); 
so that dAc = a dPc- 

Note 2. — The condition that the average velocity of flow should 
be as nearly as possible independent of the level of water in the 
canal is important for two reasons. At low velocities suspended 
matter is liable to settle on the bottom, especially in the case of a 
sewer. If the velocity is too high the material of the walls of the 
canal is liable to be washed out. Unless the shape of the cross 
section is specially proportioned, the velocity is liable to be insuffi- 
cient at low levels and too great at high levels. 

Prob. 26. — Show that the condition dAc = a dPc, deduced in 
the preceding problem also holds true for a more general form of 
exponential curve 

d=x = ^(626«'/'' + a2e->'/«), (7) 

of which the curve of Prob. 25 is a special case, when a = 6. Note. 
— The same curve expressed in hyperbolic fimctions is 

y = a [cosh-i (x/a) — cosh"^ (c/a)\ .... (8) 

where c is a constant which may be expressed through h (see Prob. 
28). The problem is more easily solved using hyperbolic fimctions. 
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Prob. 27. — Show that in Fig. 49 the lower part of the canal 
need not be a semicircle, but may be any section whatsoever, pro- 
vided that the ratio of the cross section to the wetted perimeter of 
the lower part is equal to coefficient a in Eqs. (6), (7), and (8). 
Solution. — Let Ai and Pi be the area and the wetted perimeter of. 
the lower part, and let A„ and Pu refer to the upper part. The 
condition that the average velocity be constant is 

#±4^ = const. 
Pi + P„ 

This equation must hold true at the lower limit when A« and P« are 
equal to zero. Hence 

Ai + Au __ Aj 

Pl+Pu Pi' 

or, after reduction, Ai/Pi = A^/Pw But by supposition the 
upper part of the section is shaped by a curve such that A„ = aP„ 
(see Probs. 25 and 26 above). Hence, Ai/Pi = a, which was to be 
proved. 

Prob. 28. — Deduce the most general equation for the upper 
curve in Fig. 49, to satisfy the condition dA = a dP. It has been 
shown in the preceding problems that with this condition fulfilled 
the average velocity of the water is practically independent of its 
level. Sol ution. — The g iven condition may be written in the form 
xdy = a Vl + {dx/dyY • dy, or canceling dy and squaring, 
X* = a* + a* {dx/dyy. Retaining the plus sign only, we get 

J adx /r.\ 

^2/ = -7=== (9) 

V x* — a* 

Integrating this, we obtain 



2/ = alog(x+ Vx2-a2) + const. . . . (10) 

Let the width of the cross section at LL be 2 c, and let the new 
origin be at point 0. Then the constant of integration is deter- 
mined by the condition that y = when x = c. Therefore 

Introducing, for the sake of abbreviation, 



h^c+V(?-a\ (12) 

/Google 
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and getting rid of the logarithm, Eq. (7) is obtained, which is thus 
the most general possible form of the curve. Note, — The integral 
of Eq. (9) in hyperbolic functions \a y = a cosh-^ (x/a) + const., 
the constant being determined from the condition that 2/ = when 
x = c. Eq. (8) then follows directly, thus showing again the advan- 
tage in using hyperbolic functions. 

Problems on Minimum Cost of Excavation. 

When water is used for power or irrigation purposes, 
it is sometifnes necessary to divert part of a stream at 
a smaller slope so as to create a hydraulic head or 
difference of levels (Fig. 50). Let AB represent the 
longitudinal profile of the siuf ace of a stream having a 
grade /, and let it be required to have some water power 




- A river with a derivative canal at a lesser slope. 

at C, at a head h above the stream, the rate of con- 
sumption of water or the discharge being Q. The 
water is to be obtained by means of an open canal 
branching from the stream at D and running to C at 
a lesser grade i.* The problem is to determine point 
D so that the cost of excavation of the canal, or the 
amount of material to be removed, will be a minimum. 
That there is such a minimum will be seen by tak- 
ing a point D' further up stream. Canal CD' would 
be longer than CDj but would have a smaller cross 
section since its grade is greater and consequently 
water would flow through it at a higher velocity. On 
the contrary, taking a point nearer than D, the length 
of the canal is reduced, but the slope is reduced also; 

* For a similar treatment involving closed pipes see Problems 48 
and 49 in the preceding chapter. 
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consequently the canal would have to have a larger 
cross section, and the volume of excavation may be 
larger than for a longer canal. 

Let us consider first the case of a rectangular canal 
so shaped as to possess the most advantageous cross 
section. According to Prob. 1, if the depth of such 
a canal is y, its width is 2 y, so that Eq. (3) becomes 

Q = C'2y^Vli^ (13) 

The length of the canal is determined by its slope; 
namely, the difference of the levels between the river and 
the canal increases by (7 — i) per unit length. There- 
fore 

L(I -i) =h, (14) 

where L is the length of the canal. The condition of 
minimum excavation is 

AL = ^^ = mm (15) 

I — I 

Substituting the value of y from Eq. (13) and combin- 
ing all the constants into one, we get 

const. 
-7 = nun. 

i* (7 - i) 

In order that this expression may be a minimum, its 
denominator must be a maximum, or 

i* (7 - i) = max. .... (16) 

Differentiating with respect to i and equating the 
derivative to zero, we get 

fH(7-i) -t1 = 0. 

Multiplying this equation by i* and solving for i we 
obtain 

i = f/. ..... . (17) 
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Thus, the most economical grade for the canal is equal 
to J of that of the stream. The result can also be ob- 
tained without dijBferentiation by using theorem 4 in 
the Appendix. 

Prob. 29. — Investigate the influence of the grade upon the cost 
of excavation, that is, plot to values of grade around ^ as abscissse 
the corresponding values of length, cross section, and volume of 
excavation. Assume the values corresponding to the most econom- 
ical slope to be equal to 100 per cent. 

Ana. i/I = 0.2; cost = 103.1 per cent; 
i/I = 0.4; cost = 104 per cent. 

Prob. 30. — Show that Eq. (17) is also true for a semicircular 
canal. 

Prob. 31. — Show that Eq. (17) is true for a trapezoidal section 
satisfying the condition expressed in Prob. (7). 

Prob. 32. — Show that Eq. (17) holds true for any cross section 
whatsoever, provided that it is completely determined by one of 
its dimensions. Hint, — Let z be this dimension; then A = w«*, 
P = m. 

Prob. 33. — Determine the most economical longitudinal grade 
for a rectangular canal when its depth y is given, so that the condi- 
tion X = 2y cannot be fulfilled. Hint, — It is easier in this case to 
solve Eq. (3) for i; substitute its value into the equation corre- 
sponding to Eq. (15). Having found the width x, i can be found 
from Eq. (3). 

Ana. The equation for x is /cVx* = 8 Q*y + 3 Q*x. 
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APPENDIX. 

In many practical problems the product of two 

variables must be constant, while their sum must be a 

minimum. The mathematical problem is therefore to 

find the values of u and v that satisfy the conditions 

uv = C = const.; 

u + v = min. 

Substituting the value of v from the first equation into 
the second, we obtain 

u + C/u = min. 
Equating the first derivative to zero, we get 
1 - C/u^ = 0, 

OT u = VC; consequently v = VC. We thus arrive 
at the following result : 

Theorem 1. — When the product of two variables is 
constant^ their sum is a minimum when the variables are 
equal to each other. 

Knowing this theorem, it is not necessary to differ- 
entiate in each particular problem, but the solution 
may be written directly. See, for example, the third 
solution to Prob. 1 in Chapter VI. 

Sometimes the variables satisfy a more general con- 
dition 

yMyn = (7 = const., (1) 

and the expression to be made a minimum is of the 

form 

i^p + t;« = min., (2) 

where m, n, p and q are known constants. 
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100 APPENDIX 

The total differential of the first expression is equal to 
zero, because the expression itself is a constant; the 
total differential of the second expression is also equal 
to zero, because the expression must be a minimum. 
We thus have 

Tpu^~^ du + gt;«~^ dv = 0. 
Eliminating du and dt;, and simplifying, we get 

v^ np 

The values of u and v are obtained by solving Eqs. (1) 
and (3) as simultaneous equations. This result we 
shall call 

Theorem 2. — When two variables satisfy condition 
(1), and expression (2) must be a minimum, the values of 
the variables must satisfy Eq. (3). 

In some problems the sum of two variables is con- 
stant, and their product must be a maximiun; for in- 
stance, when it is required to find a rectangle of max- 
imum area, the perimeter being given. The conditions 

are: 

u +v = const. ; 

uv = max. 

Proceeding as in the proof of Theorem 1, we get 

Theorem 3. — When the sum of two variables is con- 
stant, their product is a maximum when the variables are 
equal to each other. 
In a more general case we may have: 

t^« + yn = const., (4) 

u^v^ = max (5) 

Proceeding as in the case of Theorem 2, we derive the 
following theorem: 
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APPENDIX 101 

Theorem 4. — When two variables satisfy condition 
(4), and expression (5) mvst be a maximum , the values of 
the variables mvst be such that 

h:=m (6) 

V" mq 

Result (17) in Chapter VI may be obtained directly 
from Eq. (16) on the basis of this theorem. The 
student will find throughout the book a number of 
other problems on maxima and minima in which the 
actual differentiation can be avoided by using one of 
these four theorems. 

Strictly speaking, in order to complete the proofs 
it is necessary to show that the critical values of the 
variables found above correspond to a minimum or a 
maximum respectively. This may be done by consider- 
ing the sign of the second derivative, by interpreting 
the results geometrically or perhaps by some special 
reasoning. 
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